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Preface 


This  project  has  provided  me  with  endless  intrigue.  In  considering 
the  problem  of  guiding  a  hyper-velocity  missile  from  a  moving  platform  to 
a  moving  target,  a  common  perception  said  it  was  just  another  Kalman 
Filter  design  challenge.  I  began  with  lofty  intentions  of  running  a  full 
up  extended  Kalman  filter  simulation  but  I  became  enthralled  with  the 
limited  work  available  on  the  physical  limitations  of  basic  guidance 
algorithms.  This  became  the  emphasis  of  this  project.  The  best  filter 
scheme  cannot  overcome  a  guidance  algorithms  physical  limitation.  Only 
with  an  understanding  of  these  limitations  can  one  effectively  integrate 
a  guidance  scheme  with  a  Kalman  filter.  This  effort  is  an  exploration  of 
these  physical  limits  for  four  common  guidance  algorithms  proposed  for  the 
HVM.  My  most  untimely  PCS  to  Griff iss  AFB,  and  away 

from  the  AFIT  resources,  forced  the  construction  of  an  Extended  Kalman 
Filter  and  the  full  error  sensitivity  analysis  to  be  a  recommended  future 
effort.  For  the  investigation  and  documentation  I  completed  I  thank  my 
loving  wife,  who  has  kept  my  perspectives  straight  and  my  nose  to  the 
grind  wheel.  I  also  thank  my  God  and  Savior,  who  has  told  me  "Whatsoever 
ye  do,  do  all  to  the  glory  of  God"  [  I  Corinthians  10:31  ] 
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ABSTRACT 


This  analysis  is  an  examination  of  four  guidance  algorithms  proposed  to 
guide  a  hyper-velocity  missile  (HVM)  from  a  launch  aircraft  to  a  ground 
target.  A  technically  risky  flight  demonstration  program  of  the  HVM  on  an 
AFTI/F-16  has  been  proposed.  The  four  algorithms  are  the  Line  of  Sight 
(LOS)  (Often  called  the  Beam  Rider),  the  Line  of  Sight  Plus  (LOS+) ,  the 
Pursuit,  and  the  Proportional  guidance  algorithms.  A  simulation  of  the 
HVM  is  used  to  determine  the  no-noise  capability  of  the  algorithms  within 
the  HVM  weapon  system.  The  algorithms  are  each  developed  to  exploit  the 
launch  aircraft's  sensors  to  track  the  target  and  the  missile,  with  only 
turn  commands  being  sent  to  the  missile.  Each  algorithm  is  tested  in 
several  scenarios,  including  one  where  a  FLIR  is  the  only  tracker  used. 
A  comparative  analysis  of  the  four  algorithms  is  then  accomplished.  The 
analysis  shows  that  the  LOS  (Beam  rider)  algorithm  is  the  best  suited 
algorithm  for  the  air  to  ground  guidance  of  the  HVM.  The  proportional 
guidance  algorithm  is  an  excellent  contender  but  it  exhibits  a-  strong 
dependence  on  range  and  range  rate  measurements  from  the  launch  aircraft, 
i.e.  an  accurate  radar  must  be  used.  The  results  of  this  analysis  are 
preliminary  since  a  missile  tracking  error  analysis  has  not  been 
accomplished  to  date.  In  the  absence  of  FLIR  tracking  errors  and  system 
noise  the  LOS  guidance  algorithm  can  effectively  guide  a  HVM  to  a  direct 
hit  on  a  moving  ground  target. 
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COMPARATIVE  ANALYSIS  OF  GUIDANCE  ALGORITHMS 
FOR  THE  HYPER-VELOCITY  MISSILE  AND  AFTI/F-16 


I.  Introduction 

1 . 1  Background 

The  flight  demonstration  program  of  the  Hyper-velocity  Missile  (HVM) 
using  the  AFTI/F-16  as  the  launch  platform  has  been  proposed  and  its 
technical  risk  is  currently  assessed  as  high  [1] .  Some  major  areas  of 
concern  in  this  weapon  integration  process  are  Forward  Looking  Infrared 
(FLIR)  tracking  difficulties,  computational  loading  of  the  aircraft  fire 
control  computer,  the  HVM  performance  effects,  and  the  uncertainties  of 
missile  update  communications.  The  HVM  is  designed  to  be  a  low  cost, 
light  weight  air  to  ground  weapon.  To  achieve  this,  it  is  built  with  a 
minimal  guidance  capability  and  no  target  seeker.  Consequently,  the  HVM 
relies  on  the  launch  platform  to  guide  it  into  the  target.  The  launch 
aircraft  is  required  to  (1)  track  the  moving  target,  (2)  track  the  flying 
missile,  (3)  compute  a  new  guidance  solution,  and  (4)  communicate  this 
solution  to  the  airborne  missile.  This  situation  is  made  technically 
risky  by  the  fact  that  the  HVM  flies  so  fast  that  it  has  only  2  to  4 
seconds  of  flight  time,  and  the  very  small  missile  has  to  be  tracked  with 
a  FLIR  that  may  not  see  well  through  the  missile's  rocket  plume,  or 
through  atmospheric  conditions. 

Tht  ::omplete  magnitude  of  these  technical  problems  is  not  fully 
understood.  Efforts  have  to  be  initiated  to  eliminate  or  greatly  reduce 
these  technical  risk  areas  before  the  AFTI/F-IG  and  HVM  integration 
proceeds.  A  first  step  towards  reducing  risk  is  the  development  cf  a 
weapon  system  model  or  simulation  and  the  exploration  of  system 
sensitivities  in  these  problem  areas. 
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1 . 2  Problem  Statement 

With  several  uncertainties  in  the  HVM  weapon  system  concept  an  initial 
exploration  of  guidance  schemes  and  algorithms  is  required  to  define 
sensor  and  missile  requirements.  The  effects  of  removing  the  target 
tracker  from  the  missile  and  of  tracking  the  m.issile  and  the  target  from 
a  stand-off  flying  platform  are  not  fully  understood  or  tested.  The 
guidance  requirements  for  a  HVM  which  accepts  only  steering  commands  with 
minimal  on-board  processing  need  better  definition.  These  guidance 
algorithms  will  require  extreme  accuracy,  and  an  efficiency  that  will  not 
expend  the  missiles  limited  control  resources .  The  purpose  of  this 
research  project  is  to  develop  and  build  a  simple  weapon  system  model,  and 
perform  preliminary  sensitivity  analysis  of  the  proposed  guidance 
algorithms . 

1 . 3  Project  Scooe 

This  project  includes  the  construction  of  a  weapon  system  simulation 
and  the  analysis  of  four  HVM  guidance  techniques  against  launch  geometry 
effects  and  against  range  and  range  rate  requirements.  The  guidance 
techniques  explored  are  the  line  of  sight  (LOS)  guidance,  the  modified 
line  of  sight  (LOS+)  guidance,  the  pursuit  guidance  and  the  proportional 
guidance  algorithms.  This  selection  encompasses  the  current  guidance 
schemes  in  use  today.  This  project  models  only  one  missile  in  flight,  and 
sensor  and  system  noise  is  not  modeled  since  only  the  geometry  effects  are 
studied.  The  exploration  of  the  launch  envelope  includes  30  degree 
off-axis  shots  and  launch  during  a  5  G  aircraft  maneuver  Shots  of  12,000 
feet  range  are  tested,  from  a  F-16  model  flying  650  fps  at  1000  feet 
altitude.  It  does  not  include  maximum  or  minimum  range  shots  because  they 
are  not  certain  in  the  current  HVM  design.  A  3DOF  model  is  used  for  the 
target.  F-16,  and  HVM  so  AOA  and  sideslip  issues  are  not  addressed. 
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1 . 4  Summary  of  Current  Knowledge 

1.4.1  Introduction .  The  purpose  of  this  section  is  to  investigate  work 
related  to  modeling  the  AFTI/F-16  and  Hyper-velocity  Missile  (HVM)  weapon 
system  and  to  determine  the  current  status  of  the  HVM  to  aircraft 
integration  technology.  The  HVM  is  a  low  cost,  light  weight  anti-tank 
weapon  that  is  designed  to  increase  the  kill  ratio  of  a  modern  tactical 
air  vehicle.  In  the  attempt  to  make  the  missile  as  simple  as  possible 
most  of  the  guidance  ^responsibility  has  been  delegated  to  the  carrier 
aircraft,  requiring  some  new  technology  and  state-of-the-art  integration. 
This  review  is  an  attempt  to  determine  the  current  status  of  this 
integration  challenge  and  determine  the  feasibility  of  effectively 
modeling  the  HVM  guidance  alaorithms  to  be  implemented  by  the  AFTI/F-16. 

To  better  understand  the  integration  requirements,  some  articles 
concerning  the  HVM  and  its  current  status  are  examined  first.  Next,  two 
Master's  theses  applicable  to  this  problem  are  reviewed  as  a  first  step  in 
the  understanding  of  the  integration  issues.  This  literature  review  shows 
that  research  work  has  been  done  on  the  HVM  integration  but  there  has  been 
no  analysis  of  the  proposed  HVM  guidance  solutions  for  the  current  missile 
design . 

1.4.2  The  HVM  Concept .  The  HVM  concept  has  matured  since  T.C.  Aden 
explained  it  in  his  Hyper-velocity  Missile  paper  [2],  but  his  work  still 
details  the  heart  of  the  weapon  system.  Aden's  technical  description  of 
the  HVM  and  the  aircraft  requirements  for  an  effective  launch  are  altered 
only  because  some  of  the  missile  guidance  and  communication  details  have 
been  revised.  His  paper  presents  the  historical  problems  of  the  Close  Air 
Support  (CAS)  and  battlefield  interdiction  using  expensive  and  heavy 
conventional  precision  guided  munitions  (i.e.,  MAVERICK  and  HELLFIRE 
missiles)  and  then  presents  the  HVM  as  an  attractive,  affordable  solution 
for  multiple  target  kills  in  one  pass. 
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The  concept  of  operation  is  discussed  in  great  detail  with  excellent 
coverage  of  the  Active  Electro-Optical  Guidance  System  (AEOGS) .  The 
responsibilities  of  the  carrier  aircraft  and  its  AEOGS  are  outlined  very 
well  by  Aden  [2]  ;  however,  some  critical  emerging  technologies  are  glazed 
over  without  much  detail.  The  first  of  these  is  the  ability  of  the  AEOGS 
to  place  several  LASER  rasters  within  the  field  of  view  of  the  flying 
missiles,  a  course  raster  for  all  airborne  missiles  and  a  fine  raster, 
centered  on  the  target,  for  each  individual  missile.  (The  HVM  concept  is 
a  multi-target  kill  with  track-until-impact  requirements.)  Secondly,  the 
missile  has  to  interpret  this  raster  into  a  position  update  with  timing 
information  alone.  This  requirement  entails  some  very  fine  time 
synchronization  that  is  not  often  achieved  on  airborne  platforms,  but  will 
directly  affect  the  probabilities  of  kill  for  an  individual  missile. 

The  article  then  examines  the  HVM  in  greater  detail  with  an  emphasis  on 
the  guidance  and  control  concept.  The  guidance  and  control  concept 
explained  by  Aden  was  abandoned  during  early  development  testing  according 
to  several  conversations  with  Mr.  Boone  of  the  AFTI/F-16  SPO  and  Mr. 
Reilly  of  the  Eglin  AFB,  HVM  SPO.  Off-angle  problems  with  the  beam  rider 
technique  and  timing  problems  with  the  laser  raster  scheme  led  to  a 
current  concept  where  the  aircraft  performs  all  guidance  calculations  and 
sends  updates  to  the  missile  in  flight. 

Aden' s  paper  provides  good  coverage  of  most  of  the  HVM  concepts  but  it 
is  dated  because  of  the  new  guidance  concept  and  increased  carrier 
platform  responsibilities.  Although  this  type  of  conceptual  data  is 
required  to  construct  a  functional  model  it  is  only  a  first  step. 

1.4.3  The  HVM  Development .  The  DMS  Market  Intelligence  Report  on 
missiles  [3]  was  consulted  for  the  current  status  of  the  Hyper-velocity 
Missile  development.  DMS  Market  Intelligence  Reports  are  periodicals  that 
give  the  latest  status  of  most  DOD  research  and  development  projects.  The 
report  on  the  HVM  is  dated  Sep  1988  and  provides  a  very  compact  but 
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informative  outline  of  the  development  phases.  It  consists  of  two 
sections  entitled  Data  and  Analysis. 

The  Data  section  of  the  DMS  report  is  an  outline  of  the  pertinent 
program  data.  It  lists  the  Armament  Test  &  Development  Center,  at  Eglin 
AFB,  FL,  as  the  executive  developer  and  LTV  Corporation,  Missiles 
Division,  Dallas,  Texas,  as  the  manufacturer.  The  data  section  lists 
possible  launch  platforms,  gives  the  missile  characteristics  and  provides 
good  coverage  of  the  development  timetable  and  forecast.  It  lists  the 
ground  test  demonstration  completed  in  June  of  1987  and  a  flight  test 
demonstration  completion  in  April,  1989,  with  a  FSD  decision  in  FY  '92. 

The  Analysis  section  is  a  general  system  description  followed  by 
program  history,  plans,  current  status  and  outlook.  This  section  also 
covers  several  related  development  activities  and  gives  a  DMS  analysis  of 
the  program.  The  system  description,  although  not  extremely  technical, 
makes  reference  to  a  Multifunctional  Infrared  Coherent  Optical  Scanner 
(MICOS) ,  a  C02  laser  sensing  device.  The  MICOS  is  said  to  receive  pulses 
containing  position  information  that  must  be  time  translated  into  the 
missile  aim  point.  It  is  not  clear  from  this  report  what  changes  to  the 
guidance  scheme  were  implemented  during  ground  testing,  but  the  report 
says  the  program  has  good  joint  service  support  and  will  proceed  to  the 
full  scale  development  decision.  That  decision  will  require  information 
from  system  models. 

1.4.4  Past  HVM  Modeling.  In  1981  Captain  Donald  Capps  and  Captain 
Donald  Nelson,  in  a  MS  thesis,  developed  a  systematic  procedure  to  design 
Kalman  Filters  for  the  HVM  system  [4]  .  This  pioneer  research  work 
developed  a  workable  three  degree-of-f reedom  (DOF)  model  of  the  HVIi  and 
applied  the  latest  Kalman  Filter  modeling  tools  to  the  problem  of  guiding 
this  missile  to  a  target.  Capps  and  Nelson  assumed  a  smart  missile,  with 
a  proportional  guidance  that  did  not  account  for  the  aircraft  calculating 
a  missile  to  target  line-of-sight  rate.  Although  this  research  speculated 
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about  some  of  the  missile  parameters  and  the  missile  guidance  technique, 
it  is  a  major  asset  in  developing  much  of  the  weapon  system  model. 

Their  work  put  much  effort  in  building  an  accurate  truth  model  of  the 
HVM.  It  also  evaluated  some  line-of-sight  and  inertial  navigation  system 
Kalman  Filters  in  an  attempt  to  come  up  with  a  reduced  order  Kalman  Filter 
that  could  accurately  estimate  the  system  performance.  A  line  of  sight 
Kalman  filter  was  determined  as  the  best  method  of  tracking  the  HVM 
missile.  This  is  the  tracking  method  to  be  used  in  this  guidance 
algorithm  analysis. 

The  work  done  in  Capps'  thesis  used  guidance  concepts  that  required  a 
guidance  system  within  the  missile,  and  assumed  a  proportional  guidance 
algorithm.  Although  it  is  a  major  stepping  off  point  for  the  full 
development  of  an  AFTI/F-16  and  HVM  weapon  system  model,  it  did  not 
e.xplore  the  guidance  solutions  applicable  for  the  current  HVM 
developments . 

1.4.5  AFTI/F-16  Modeling  Effort.  Captain  David  Michalk,  in  a  MS  thesis, 
investigated  the  target  state  estimation  techniques  for  the  AFTI/F-16 
automated  maneuvering  and  attack  system  [5]  .  This  investigative  work 
provides  excellent  coverage  of  the  AFTI/F-16  sensor  suite,  the  target 
state  estimator,  and  the  on-board  implementation  of  a  Kalman  Filter. 
Although  this  work  was  done  for  the  Air-To-Air  mode,  the  details  of  the 
document  provide  many  of  the  current  modeling  techniques,  and  an  insight 
to  unique  AFTI/F-16  capabilities  that  can  directly  relate  to  the 
Air-to-Ground  mode  modeling  for  an  HVM  analysis.  It  is  expected  that  the 
AFTI/F-16  FLIR  capabilities  will  be  strained  attempting  to  track  a  small 
fast  missile  that  is  in  the  line  of  sight  of  the  target.  No  other 
documentation  addressing  this  problem  appears  to  be  available. 

1.4.6  Review  Conclusions.  Much  work  has  been  accomplished  in  the 
modeling  of  AFTI/F-16  and  HVM  systems.  It  is  necessary  to  pull  this  work 
together  into  a  weapon  system  model  that  explores  pertinent  guidance 
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algorithm  sensitivities  before  base  line  fire  control  algorithm.s  are 
developed.  Research  was  conducted  during  the  conceptual  phase  of  the  HVI-^ 
development,  but  there  is  no  current  analysis  that  fully  explores  the 
risky  integration  issues.  The  HVM  flight  demonstration  program  may  soon 
be  underway,  and  the  AFTI/F-16  to  HVM  integration  effort  needs  some  up¬ 
front  research  to  select  a  guidance  scheme.  There  are  related  analysis 
efforts  with  spin-offs  that  might  be  directly  applicable  to  such  an  effort 
and  the  modeling  tools  required  to  perform  this  analysis  are  in  place. 
Further  research  with  the  latest  missile  parameters  and  the  current 
guidance  techniques  can  greatly  enhance  this  integration  process  and  more 
fully  demonstrate  the  effectiveness  of  this  weapon  system. 

1 . 5  Approach  Methodology 

1.5.1  Introduction .  This  research  takes  the  first  step  of  an 
algorithmic  analysis.  The  examination  of  the  algorithms  and  the  analysis 
accomplished  amount  to  the  preliminary  evaluation  of  the  candidate 
algorithms  under  no-noise  conditions.  This  step  is  important  because  a 
guidance  algorithm  can  do  no  better  than  the  controlling  physics  involved. 
Even  at  this  level  the  differences  discovered  in  the  guidance  algorithms 
performance  for  a  hyper-velocity  missile  are  striking  and  informative.  A 
more  in-depth  analysis  with  both  system  noise  and  measurement  noise  is 
required  on  the  algorithm  selected  after  this  cursory  analysis. 

The  approach  taken  in  the  development  and  analysis  of  the  weapon  system 
guidance  algorithms  may  have  some  unique  situations  and  parameters  but  it 
is  a  standard  weapon  system  modeling  approach.  It  consists  of  building  an 
aircraft,  missile,  and  target  model  that  interact  dynamically  to  provide 
a  weapon  system  simulation.  The  simulation  begins  with  the  use  cf 
position  truth  models  for  the  launch  aircraft,  the  target  vehicle  and  the 
missile.  A  radar  and  FLIR  tracker  model  takes  this  true  position  data  and 
calculates  the  line  of  sight  and  range  parameters  that  would  be  available 
to  the  guidance  algorithms.  The  selected  guidance  algorithm  then  provides 
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the  guidance  command  data  to  the  missile  truth  model.  The  missile  truth 
model  applies  this  command  to  a  second  order  response  and  continues  the 
missile  flight  until  ground  impact.  The  miss  distance  to  the  target  is 
the  measure  of  effectiveness  for  each  simulation.  A  block  diagram  of  this 
simulation  is  shown  in  Figure  1-1.  It  is  explained  in  greater  detail 
below  with  a  discussion  of  the  analysis  effort  undertaken  in  this 
research . 

1.5.2  The  Truth  Model.  The  truth  model  calculates  the  position  and 
velocity  of  the  aircraft,  the  target,  and  the  missile.  The  aircraft  model 
uses  the  scenario  setup  variables  to  establish  a  three  dimensional,  three- 
degree-of-f reedom  (3D0F)  flight  profile  for  the  aircraft.  The  target 
model  does  the  same  for  the  target  but  its  vertical  plane  is  held 
constant.  The  missile  model  is  of  higher  fidelity  than  the  aircraft  or 
target  models.  The  missile  acceleration  calculations  must  include  the 
rocket  motor  thrust,  air  drag,  lateral  acceleration  commands,  and 
aerodynamic  turning  drag.  Once  these  are  calculated  the  missile  model  also 
computes  the  current  missile  position  and  velocity  in  the  three 
dimensional,  three-degree-of-f reedom  reference  system. 

The  truth  models  take  scenario  inputs  and  missile  parameters  from  the 
Scenario  Setup,  and  missile  commands  are  received  from  the  selected 
guidance  algorithm.  Aircraft,  target,  and  missile  position  and  velocity 
are  sent  to  the  measurement  model  at  regular  time  intervals  established  by 
the  Scenario  Setup. 

1.5.3  The  Measurement  Model.  The  measurement  model  implemented  in  this 
analysis  is  a  simple  coordinate  transformation  that  calculates  a  relative 
LOS,  LOS  rate.  Range  and  Range  rate  from  the  aircraft  to  both  the  target 
and  the  missile.  The  LOS  and  LOS  rate  are  calculated  in  the  horizcr.tal 
plane  as  azimuth,  and  in  the  vertical  plane  as  pitch.  The  measurerr.ent 
model  provides  identical  data  to  each  of  the  algorithms  under  test  with  a 
flavor  of  realism.  For  this  initial  and  simple  analysis,  signal  noise  is 
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not  accounted  for,  and  the  construction  of  a  Kalman  filter  that  would 
smooth  the  actual  measurement  data  is  not  accomplished.  The  accurately 
calculated  truth  data  is  initially  input  into  the  guidance  algorithm  to 
establish  the  upper  limit  of  performance. 

The  measurement  model  uses  the  inertial  position  and  velocity  data  of 
the  aircraft,  the  target,  and  the  missile.  It  then  outputs  the  relative 
LOS,  LOS  rate,  Range,  and  Range  rate  data  to  the  selected  guidance 
algorithm.  The  time  rate  at  which  the  measurement  model  performs  the 
transformation  is  established  by  the  Scenario  Setup. 

1.5.4  The  Guidance  Algorithm.  The  guidance  algorithm  takes  the  output 
from  the  measurement  model  and  uses  this  data  to  calculate  the  missile 
commands  sent  to  the  missile  truth  model.  A  guidance  algorithm  must  be 
chosen  at  initiation.  LOS  guidance,  LOS+  guidance.  Pursuit  guidance  and 
Proportional  guidance  methods  are  available.  These  guidance  methods  are 
discussed  in  detail  in  Chapter  3- 
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1 . 6  Surronarv 


This  analysis  is  an  examination  of  four  guidance  algorithms  proposed  to 
guide  a  hyper-velocity  missile  (HVM)  from  a  launch  aircraft  to  a  ground 
target.  A  technically  risky  flight  demonstration  program  of  the  HVM  on  an 
AFTI/F-16  has  been  proposed.  Chapter  2  describes  the  system  models 
developed  to  test  the  guidance  schemes.  Chapter  3  develops  the  Line  of 
Sight  (LOS)  (Often  called  the  Beam  Rider),  the  Line  of  Sight  Plus  (LOS+) , 
the  Pursuit,  and  the  Proportional  guidance  algorithms.  A  simulation  of 
the  HVM  is  used  to  determine  the  no-noise  capability  of  the  algorithms 
within  the  HVM  weapon  system.  The  algorithms  are  each  developed  to 
exploit  the  launch  aircraft's  sensors  to  track  the  target  and  the  missile, 
with  only  turn  commands  being  sent  to  the  missile.  Chapter  4  documents 
the  testing  of  each  algorithm  in  several  scenarios,  including  one  where  a 
FLIR  is  the  only  tracker  used.  A  comparative  analysis  of  the  four 
algorithms  is  then  accomplished.  The  conclusions  and  recommendations  are 
shown  in  Chapter  5. 

With  an  understanding  of  the  background,  the  scope  and  the  approach 
methodology  it  is  important  to  understand  the  system  truth  models 
developed  in  the  next  chapter. 
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II.  SYSTEM  TRUTH  MODEL 


2 . 1  Introduction 

The  computer  simulation  of  a  system  is  always  a  compromise  of 
fidelity.  In  the  case  of  a  weapons  system  model  it  is  important  that  the 
truth  model  have  greater  fidelity  than  the  algorithms  under  test.  The 
model  used  for  this  simulation  contains  a  position  and  velocity  truth 
model  that  uses  a  double  precision  Fortran  routine.  The  other  portions  of 
code,  those  that  might  be  implemented  in  an  airborne  computer,  are 
implemented  in  single  precision  Fortran  code.  To  attain  greater  precision 
in  the  truth  model  the  integration  time  is  kept  much  smaller  than  the 
sensor  or  the  missile  time  responses.  Figure  1-1  is  a  block  diagram  of 
the  general  system  that  is  built  and  used  for  preliminary  tests  of  the 
four  proposed  algorithms.  Figure  2-1  shows  the  geometric  layout  of  the 
model  scenarios.  Since  the  simulation  uses  a  three  dimensional  space 
there  is  a  synonymous  geometry  in  the  vertical  plane,  which  is  shown  in 
Figure  2-2.  The  coordinate  system  used  in  this  work  is  the  geographic 
North,  East,  Down  (NED)  coordinate  system  as  seen  in  these  Figures.  The 
order  of  sections  in  this  chapter  are  the  assumptions  and  detailed 
missile,  aircraft,  and  target  truth  models,  and  the  calculations  necessary 
to  construct  a  measurement  model.  The  purpose  of  constructing  these  truth 
and  measurement  models  is  to  analyze  the  guidance  algorithms  described  in 
the  next  chapter. 

2 . 2  Model  Executive  and  Assumptions 

This  model  is  assembled  in  as  simple  a  manner  as  possible  while 
affording  a  versatile  analysis  capability.  This  noise  free  model  is 
initially  constructed  as  a  tool  to  test  the  guidance  algorithm  code. 
After  only  a  little  testing  of  the  algorithms  it  became  clear  that  a  clean 
and  simple  weapon  system  model  is  adequate  for  the  very  preliminary 
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Figure  2-1  Weapon  Model  Horizontal  Plane  Geometry 


Figure  2-2  Weapon  Model  Vertical  Plane  Geometry 


algorithm  analysis  that  has  yet  to  be  accomplished  for  a  hyper-velocity 
missile,  monitored  and  controlled  from  the  launch  platform.  Thus  the 
model  is  a  simple  but  accurate  implementation  that  enables  most  of  the 
research  effort  to  be  centered  around  the  guidance  algorithms. 

The  Scenario  Setup  block  shown  in  Figure  1-1  allows  versatile  control 
of  a  missile  launch  scenario.  The  use  of  the  I/O  subroutine  allows 
missile  variables  and  scenario  variables  to  be  adjusted  from  default 
values  before  each  run.  The  system  timing  and  the  guidance  algorithm  is 
selected  and  controlled  via  the  scenario  variable  selection. 

The  truth  models  are  stepped  through  their  motion  at  a  small  time 
interval.  (The  default  time  delta  is  set  to  2  milliseconds.)  On  a  single 
pass  the  system  time  is  incremented  by  DT  (the  system  time  interval)  ; 
missile,  aircraft,  and  target  truth  models  are  updated;  a  missile 
proximity  to  target  check  is  made;  excessive  flight  time  and  missile  range 
checks  are  made;  then  the  timed  events  are  checked  for  their  execution 
time.  The  timed  events  include  the  sensor  output  updates,  the  missile 
updates  including  guidance  algorithm  execution,  and  the  recording  of 
pertinent  data.  The  timing  of  the  timed  events  is  controlled  by  scenario 
inputs.  The  missile  and  scenario  inputs  with  their  default  values  are 
shown  in  Table  2-1  and  Table  2-2  respectively. 

Throughout  the  simulation  pertinent  data  is  stored  in  a  matrix  for 
future  output  to  a  file.  Each  matrix  can  hold  ten  variables  for  200 
individual  time  samples.  This  format  is  ideal  for  future  data  analysis 
via  spreadsheets.  The  four  matrices  store  selected  control  data,  the 
missile,  aircraft  and  target  position  data,  the  missile,  aircraft  and 
target  velocity  data,  and  ten  of  the  twelve  measurement  model  states.  At 
the  end  of  a  run  the  user  is  given  the  option  of  saving  this  data  as  a 
file  or  discarding  it. 
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Table  2-1  Missile  Input  Variables 


NAME 

DEFAULT  VALUE 

UNIT 

SYMBOL 

MSL  Thrust  Time 

1 . 5 

SEC 

TBURN 

MSL  Thrust 

6000.0 

LBS 

THRUST 

MSL  Body  Weight 

25.9 

LBS 

WTBD 

MSL  Fuel  Weight 

32.17 

LBS 

WTFU 

MSL  Diameter 

3.8 

IN 

MDIA 

MSL  Length 

7.0 

FT 

MLNGTH 

MSL  Wn 

44.0 

EPS 

WN 

MSL  Damping 

0.7071 

NU 

ZETA 

MSL  Max  G 

100.0 

G 

GMAX 

MSL  Max  Fit  Time 

4.0 

SEC 

TFLGHT 

MSL  Drag  Constant 

0.000174 

NU 

KDRAG 

Table  2-2 

Scenario  Input 

Variables 

NAME  DEFAULT  VALUE 

UNIT 

SYMBOL 

Aircraft  Speed 

650.0 

FPS 

VAC 

Acft  Heading 

60.0 

DEG 

HDGA 

Acft  Altitude 

1000.0 

FT 

AALT 

Acft  Lateral  Accel 

00.0 

G 

AACf 

Acft  Climb  Accel 

0.0 

G 

APITCH 

Range  To  Target 

12000.0 

FT 

RT 

Target  Speed 

50.0 

MPK 

VTGT 

Target  Heading 

90.0 

DEG 

HDGT 

Target  Lateral  Accel  00.0 

FPS''2 

ATGT 

Sys  Integr  Time 

0.002 

SEC 

DT 

MSL  Update  Times 

0.02 

SEC 

DTMU 

Sensor  Update  Times 

0.02 

SEC 

DTSU 

Output  Times 

0.02 

SEC 

DTOUT 

Filter  Update  Times 

0.02 

SEC 

DTKU 

Correction  Time 

0 . 1 

SEC 

TC 

Other  Constant 

1.5 

SEC 

K2 

Porportional  Constant  0.003 

NU 

K3 

Algorithm  Selection 

1.0 

NU 

ALG 
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2 . 3  Missile  Model 


The  model  of  the  missile  is  an  integral  part  of  the  system  executive.  It 
consists  of  an  integration  of  the  state  vectors  at  an  integration  time  of 
DT,  the  system  time  delta.  A  block  diagram  of  this  model  is  shown  in 
Figure  2-3. 

Several  of  the  missile  parameters  are  read  from  a  data  file  called 
"msl.dat".  This  allows  for  easy  change  of  some  basic  parameters.  These 
variable  and  their  default  values  are  in  a  data  file  as  shown  in  Table 
2-1.  At  program  initiation  any  of  these  variables  may  be  altered  before 
the  program  proceeds . 

Before  the  missile  states  are  updated  in  each  cycle,  the  missile 
response  to  the  turn  command  rates  are  modeled  as  a  second  order  lag 
response.  This  second  order  model  is  controlled  by  the  missile  natural 
frequency  and  damping  coefficient  W„  and  5  respectively.  These  are  input 
from  the  data  table  and  may  be  changed  before  a  simulation.  The  default 
values  used  were  Wn=  44  radians/second  and  8  =  .707,  these  were  used  in  the 
previous  HVM  model  by  Capps  and  Nelson  [4:D-2].  The  commanded  turn  from 
the  guidance  algorithm  is  the  input  to  this  second  order  lag  and  the 
actual  turn  rate,  is  the  filter  output.  Both  the  heading  rate,  ©i^, 

and  the  pitch  rate,  0,„,  are  computed  in  the  geographic  NED  coordinate 
frame.  The  equations  for  each  iteration  are  shown  below: 


0,,  =  ©,^  + 

-  H,,,))  T 

(2-1) 

+ 

II 

(Wp," 

-  Pu,))  ^ 

(2-2) 

and 

<©in, 

- 

2 

5  Wp  H,,,)  T 

(2-3) 

Plag  + 

(4. 

- 

2 

5  Wp  P,„)  T 

(2-4) 

Where:  0^,  =  The  Actual  Turn  Rate  (radians  /  sec) 


0,^  =  The  Actual  Pitch  Rate  (radians  /  sec) 

W,,  =  Missile  Natural  Frequency  (radians  /  sec) 
6  =  Natural  Damiping  of  the  Missile  (No  Units) 
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flccel  Vel  Position 


Figure  2-3  Missile  Truth  Model  Block  Diagram 

=  The  Commanded  Turn  -ate  (radians  /  sec) 

®2cnd  “  The  Commanded  Pitch  Rate  <radiar«  /  sec) 

Hj,3,  =  Internal  2nd  order  state  variable 
Pjj,  =  Internal  2nd  order  state  variable 
T  =  System  integration  time  time  delta  (sec) 

The  missile  and  aircraft  flight  starts  at  time  zero,  location 
(x,y,z)  =  (0,0,0)  .  The  missile  has  an  initial  velocity  equal  to  the 

aircraft  velocity  and  an  acceleration  of  THRUST/MASS.  The  missile  weight 
is  decremented  during  each  update  during  the  roc)cet  burn.  The  amount  of 
the  weight  decrement  is  calculated  so  the  fuel  weight  is  eliminated  at 
burnout  time,  assuming  constant  thrust.  A  linear  fuel  burn  and  a  constant 
thrust  are  normal  missile  design  goals.  The  heading  of  the  missile  is 
altered  by  the  applied  lateral  commands.  This  is  accomplished  by  turning 
the  missile  half  the  distance  it  will  turn  in  a  time  interval  of  T  ; 
applying  the  acceleration  caused  by  the  engine  thrust;  updating  the 
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velocity  and  position  vectors;  then  turning  the  missile  the  remaining 
commanded  rotation  and  adjusting  the  velocity  vector  to  align  to  this  new 
heading.  This  method  ensures  a  smooth  and  accurate  implementation  of  the 
commanded  turn.  The  Euler  integration  technique  is  used  in  this  model  so 
it  is  important  to  keep  the  integration  step  size  much  smaller  than  the 
time  constant  associated  with  the  missile  (4:B-3].  In  this  test  a  2  msec 
integration  time  was  used. 

The  drag  on  the  missile  is  calculated  much  the  same  as  Capps  did  except 
for  a  lift  drag  penalty  proportional  to  the  turn  command  [4:2-9  -  2-10]  . 
The  drag,  n,  is  calculated  as  follows: 

H  =  -  c/2  S  (V,/V,  +  0,,  V„  K,,)/  M  (2-5) 

Where:  (1  =  Air  Drag  Estimate  (ft/sec^) 

c  =  Air  density  (.0023  slugs/ft^  ) 

S  =  MSL  Cross-sectional  Area  (ft^) 

V„  =  MSL  Velocity  (EPS) 

V3  =  Speed  of  Sound,  1087.1  (EPS) 

(i.e.  V,  /  V„  =  1/MACH#) 

0i„  =  Actual  Missile  Turn  Rate  (Radians  /  sec) 

=  A  Lift  Drag  Penalty  Constant  (sec^/ft) 

M  =  Missile  Mass  (Slugs) 

In  many  models  the  missile  drag  is  based  on  only  the  body  drag  with  no 
penalty  for  the  lift  drag  caused  by  missile  turns.  The  constant  has 
a  default  setting  which  doubles  the  air  drag  during  a  50  G  ttri  at  5000 
EPS  missile  velocity.  This  drag  penalty  is  caused  by  the  missiles  angle 
of  attack  during  turning  and  is  added  here  to  add  fidelity  to  the  HVM 
motion  model.  The  coupling  of  this  lift  drag  to  the  body  drag  through  a 
proportion  constant  is  better  understood  by  examining  BlackeJ.ock' s 
development  of  the  drag  equation  [8-332].  Thus,  in  this  model  there  is  a 
drag  penalty  added  for  high  G  turning  commands. 
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In  this  manner  the  missile  state  vector  is  updated  each  cycle  of  the 
executive.  The  state  vector  includes  the  three  dimensional  position  and 
velocity  vectors  tracked  in  double  precision.  The  calculation  of  these 
states  is  straight  forward. 

2 . 4  Aircraft  Model 

Since  the  missile  flight  is  so  short  (approximately  4  seconds  maximum) 
a  very  simple  3DOF  model  is  used  to  simulate  the  aircraft  position. 
Because  of  its  simplicity  it  is  not  a  validated  F-16  model,  but  it  can 
simulate  the  position  of  a  maneuvering  F-16  during  a  4  second  portion  of 
flight.  The  3D0F  model  of  the  aircraft  is  an  integral  part  of  the  system 
executive.  It  consists  of  an  integration  of  the  state  vectors  at  an 
integration  time  of  T,  the  system  time  delta.  A  block  diagram  of  this 
model  is  the  same  as  the  missile  model  shown  in  Figure  2-3,  except  the 
input  thrust  and  drag  are  always  equal  and  opposite.  This  assumption 
makes  the  aircraft  maintain  a  constant  speed  during  the  missile  flight, 
even  if  the  aircraft  engages  in  a  turn.  This  is  a  realistic  flight 
maneuver  for  the  short  engagement  time. 

The  aircraft  track  is  determined  by  the  variable  selection  prior  to 
the  simulation  run.  The  variables  that  may  be  altered  prior  to  each  run 
are  shown  in  Table  2-2.  The  ones  which  control  the  aircraft  model  include 
an  initial  aircraft  heading,  speed,  pitch,  turn  G  and  climb  G.  For  a  left 
turn  or  a  dive  the  turn  G  and  climb  G  are  negative.  These  initial  G 
maneuvers  are  transformed  into  a  total  G  maneuver  at  a  body  angle  to  keep 
the  coordinate  frame  proper  during  a  pitch  over  mianeuver.  Since  AOA  and 
sideslip  are  not  modeled  in  a  3D0F  model  it  is  left  to  the  care  of  the 
user  to  ensure  that  aircraft  heading  or  pitch  does  not  inhibit  sensor 
tracking  bounds  for  the  target  or  the  missile.  All  maneuvers  are  assumed 
to  occur  at  a  constant  speed,  and  the  maneuver  is  constant  throughout  the 
short  missile  flight  time.  There  is  also  no  G  onset  rate  so  when  set  for 


2-8 


a  5  G  turn  (for  example) ,  the  aircraft  is  already  executing  this  turn  at 
t=0,  launch  time. 

At  launch  time,  time  zero,  the  aircraft  position  establishes  the 
origin  of  the  reference  inertial  plane.  The  aircraft  initial  velocity 
vector  is  established  by  the  input  heading,  pitch  and  aircraft  speed.  The 
heading  of  the  aircraft  is  altered  by  the  applied  maneuver  commands .  This 
is  accomplished  by  turning  the  aircraft  half  the  distance  it  will  turn  in 
a  time  interval  of  x;  updating  the  velocity  and  position  vectors;  then 
turning  the  aircraft  the  remaining  commanded  rotation  and  adjusting  the 
velocity  vector  to  align  to  this  new  heading.  This  method  ensures  a 
smooth  and  accurate  implementation  of  the  commanded  turn,  as  discussed  in 
the  missile  model  description. 

In  this  manner  the  aircraft  state  vector  is  updated  with  the  Euler 
integration  technique  during  each  cycle  of  the  executive.  The  state 
vector  includes  the  three  dimensional  position  and  velocity  vectors 
tracked  in  double  precision.  The  calculation  of  these  states  is  straight 
forward. 

2 . 5  Target  Model 

The  target  model  used  in  this  analysis  is  a  simple  3DOF  approximation 
of  a  moving  ground  target.  The  model  employes  a  simple  Euler  integration 
technique  in  the  NED  coordinate  frame  with  an  x  and  y  (north  and  east) 
plane  and  a  constant  vertical  plane  assumed.  The  form  of  the  Euler 
integration  is  as  follows: 


V.(t+1)  =  V,(t) 

+ 

A,  dt 

(2-6) 

Vy(t+l)  =  Vy(t) 

+ 

Ay  dt 

(2-7) 

Xt(t  +  1)  =  A,  (t) 

+ 

V,  (t  +  l) 

dt 

(2-8) 

Y,  (t+i)  =  y,  (t) 

+ 

Vy(t+l) 

dt 

(2-9) 

The  initial  conditions  and  the  acceleration  terms  are  calculated  from 
the  scenario  input  variables  LOS^,  R^,  and  A^,.  shown  in  Table  2-2. 
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These  calculations  are  shown  below. 


A,  =  cos(Ht,J 
Vy(0)  =  sin(H,^) 
Yt(0)  =  Rt  sin  (LOS.) 


Ax  =  A^gt  sin(Ht,J 
Vx(0)  =  cos(H,^) 

Xt(0)  =  cos(LOSj 
Where:  LOS^ =  Line  of  Sight  angle  from  Aircraft  to  Target 
Rt  =  Range  from  Aircraft  to  Target 
=  Velocity  of  Target 
Ht,t  =  Heading  of  Target 
A^,t  =  Acceleration  of  Target 

In  this  manner  the  target  is  given  motion  for  the  algorithm  evaluation. 
The  target  integration  is  accomplished  each  cycle  of  the  execution  as  is 
the  aircraft  and  missile  Euler  integrations.  The  measurement  model  is 
used  to  calculate  the  missile  and  target  position  relative  to  the  aircraft 
and  to  put  these  in  an  angle  and  range  format  as  would  be  done  by  a  FLIR 
and  a  radar.  Since  the  aircraft  attitude  is  not  calculated  in  the  3D0F 
model,  these  lines  of  sight  measurements  are  not  in  the  aircraft  trac)cing 
frame  as  done  by  some  models.  In  the  NED  coordinate  frame  the  line  of 
sight  measurements  are  the  angle  from  north,  and  the  angle  above  the 
horizon,  as  shown  in  Figure  2-1  and  2-2. 

2 . 6  Me  lurement  Model 

The  Measurement  model  is  made  up  of  a  14  state  matrix  and  the 
calculations  required  to  update  the  matrix  from  the  inertial  data.  The 
states  selected  for  the  model  are  as  follows: 


x( 

1)  =  B,, 

=  Azimuth 

LOS 

aircraft  to  target 

x( 

2 )  = 

=  Azimuth 

LOS 

rate  aircraft  to  target 

x( 

3)  =  B,, 

=  Pitch 

LOS 

aircraft  to  target 

x( 

4)  =  Bj, 

=  Pitch 

LOS 

rate  aircraft  to  target 

x( 

It 

=  Range 

aircraft  to  target 
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x(  6)  -  R, 
x(  7)  = 
x(  8)  =  Bi„ 
x(  9)  = 
x(10)  =  Bin. 
x(ll)  =  R„ 
x(12)  = 
x(13)  = 
x(14)  = 


=  Range  rate  aircraft  to  target 
=  Azimuth  LOS  aircraft  to  missile 
=  Azimuth  LOS  rate  aircraft  to  missile 
=  Pitch  LOS  aircraft  to  missile 
=  Pitch  LOS  rate  aircraft  to  missile 
=  Range  aircraft  to  missile 
=  Range  rate  aircraft  to  missile 
=  Range  estimate,  aircraft  to  missile 
=  speed  estimate  of  the  missile  relative  to  aircraft 


These  last  two  states  are  constructed  from  a  very  simple  missile  model 
in  order  to  test  the  algorithm  with  FLIR  inputs  only,  i.e.  no  radar  range 
or  range  rate  data  available.  The  simplified  missile  model  assumes  that 
the  missile  has  constant  acceleration  and  constant  drag,  and  that  it 
flies  directly  away  from  the  launching  aircraft.  As  expected  it 
accumulates  error  rapidly,  and  this  error  is  used  to  demonstrate  guidance 
algorithm  sensitivity  to  poor  range  data. 

The  aircraft  could  use  a  Kalman  filter  to  provide  a  best  estimate  of 
these  values  to  the  guidance  algorithm.  In  this  simplified  model  single 
precision  'perfect'  measurements  approximate  the  data  that  would  be 
available  to  the  guidance  algorithms  on  the  aircraft.  The  single 
precision  used  is  an  8  bit  exponent  and  a  24  bit  mantissa.  This  gives  a 
range  accuracy  of  7x10'^  feet,  range  rate  of  4x10"^  fps,  and  angle 
measurements  of  6x10"’  rad.  Though  noise  free,  the  modeled  state 
measurements  are  adequate  for  the  this  first  level  evaluation  of  the 
guidance  algorithms. 

The  calculations  of  these  states  is  accomplished  at  an  update  rate  as 
specified  in  the  scenario  input.  The  states  are  shown  in  Figure  2-1  and 
2-2  and  are  calculated  from  the  geographic  NED  x,  y,  and  z  values  as 
follows : 
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— 

tan-M(Yt-YJ/(X,-X.)) 

(2-10) 

= 

tan'M  (Y„-yj  /  (X„-XJ  ) 

(2-11) 

rt 

= 

((X,-X.)^  +  (Y,-YJ^)'^' 

(2-12) 

Bit 

= 

tan-‘  (  (Zt-ZJ  /rj 

(2-13) 

= 

((X„-XJ^  +  (Y„-YJ^)‘'^ 

(2-14) 

B2„ 

= 

tan'M  (Z„-ZJ /rj 

(2-15) 

Rt 

= 

((Z,-ZJ"  +  r, 

(2-16) 

K 

= 

((Z„-ZJ^  + 

(2-17) 

Bu 

= 

(Vyt  -VyJ  cos(aiJ  -  (V,,  -V„  )  sin(B,J 

rt 

(2-18) 

Bi„ 

= 

(V^-V,J  cos(B:J  -  (v„-v„)  sin(BiJ 

(2-19) 

Bjt 

= 

(V,,-V„)/R,  cosdijJ  - 

(V,,-V„)/R,  cos  (a,,)  sin<Bj,)  - 

(Vy,  -Vy,) /R^  sin  (Bit)  sin(Bjt) 

(2-20) 

B,„ 

s 

<Vt„-Vt.)/R„  cos(BjJ  - 

(V„.-V..)/R„  cos(BiJ  sin(B,J  - 

(Vyf-VyJ /R„  sin(Bi„)  sin(Bj„) 

(2-21) 

Rt 

= 

(Vjt-VjJ  sin(B2t)  +  (V;,t-V^a)  cos  (Bit)  cos(B2t)  + 

(Vy.-Vy)  sin  (Bit)  cos(B2t) 

(2-22) 

K 

= 

(V„  -V„)  sin(B2J  + 

cos(B,J  cos(Bj„)  + 

(V^-VyJ  sindi.J  cos(Ii,J  (2-23) 

For  the  3D0F  implementation  of  this  model  the  missile  attitude  is  not 
calculated.  Instead,  the  missile  NED  x,  y,  and  z  values  are  taken  from 
the  missile  truth  model  and  used  to  calculate  the  parameters  needed  by  the 
guidance  algorithms . 
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2 . 7  System  Truth  Model  Summary 


The  3D0F  truth  model  constructed  for  this  analysis  is  adequate  for 
algorithm  familiarity  and  preliminary  test.  The  model  is  not  intended  for 
the  full  error  sensitiyity  analysis  where  a  6D0F  model  is  needed  for  angle 
measurement  analysis.  This  sensitivity  analysis  needs  to  be  accomplished 
on  guidance  algorithms  which,  in  this  preliminary  test,  are  most 
effective . 

The  measurement  model  is  made  to  use  the  range  and  velocity  estimates 
for  its  range  and  range  rate  measurements  during  some  tests.  This 
simulated  the  elimination  of  the  radar  sensor  and  tested  the  guidance 
algorithms  dependence  on  range  and  range-rate  measurements .  For  more 
details  on  this  analysis  see  Chapter  4. 

The  models  constructed  here  are  simple  but  useful  for  performing  the 
initial  comparative  analysis  of  the  four  guidance  algorithms.  The  next 
chapter  gives  a  detailed  explanation  of  the  guidance  algorithm 
development . 
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III.  Guidance  Algorithm  Development 


3 . 1  Introduction 

The  guidance  algorithms  tested  in  this  study  have  been  tested  in  many 
applications.  The  Line-Of-Sight  (LOS)  algorithm,  sometimes  called  a  beam 
rider  algorithm,  is  very  common  in  surface  to  air  missile  applications 
[6:158].  The  LOS+  algorithm  is  similar  to  the  LOS  algorithm  with  a 
predictor  that  estimates  the  final  line  of  sight  to  the  target  and  guides 
the  missile  down  this  line  to  the  target.  The  pursuit  algorithm  is  an  old 
standby  used  because  of  its  very  simple  implementation.  This  algorith-m  is 
implemented  by  causing  the  missile  to  continually  steer  towards  the 
target.  The  proportional  guidance  algorithm  is  a  low  control  energy 
guidance  algorithm  based  on  the  fact  that  if  the  LOS  between  a  missile  (or 
originally  a  ship)  and  a  target  does  not  change  they  are  on  a  collision 
course.  (This  is  preferred  by  the  missile-target  scenario,  but  dreaded  by 
the  ship-to-ship  scenario.) 

A  different  application  of  these  algorithms  is  represented  here  because 
they  are  being  used  in  a  system  quite  foreign  to  their  original 
development.  The  LOS  algorithm  has  been  developed  and  used  largely  for  a 
stable  platform  tracking  a  moving  target,  here  we  have  a  moving  platform 
tracking  a  relatively  stable  target.  The  pursuit  and  proportional 
guidance  algorithms  were  built  for  a  missile  sensor  tracking  a  target, 
here  a  third  party  is  involved,  an  aircraft,  trying  to  sense  the  targets 
location  relative  to  the  missile's.  In  each  application  the  guidance 
algorithms  must  deal  with  an  extremely  fast  missile  that  uses  only  turn  or 
climb  commands.  In  this  sense  then,  this  is  a  new  test  of  these 
algorithms . 

The  development  of  these  algorithms  is  discussed  in  the  sections  that 
follow.  It  should  be  remembered  again  that  the  geometry  of  the  scenario 
takes  a  little  effort  to  understand  in  each  discussion.  These  changes  in 
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the  basic  scenario  geometry  cause  the  algorithms  to  react  differently  than 
in  their  normal  implementation,  as  seen  in  the  analysis  section  of  this 
report . 

3 . 2  Line  of  Sight  Algorithm  Development 

The  LOS  guidance  principle  is  often  called  the  beam  rider  method  and 
was  developed  fo,,  ground  to  air  applications  [6:158  ] .  It  involves  keeping 
the  missile  on  a  line  projected  from  the  shooter  to  the  target.  Note  that 
this  is  different  than  a  pursuit  guidance  scheme  where  the  missile  is 
pointed  down  to  its  own  line  of  sight  to  the  target.  The  geometry  of  the 
LOS  guidance  scheme  is  shown  in  Figure  3-1. 

In  both  snapshots  of  this  figure  a  missile  steering  correction  is 
required  to  move  the  missile  onto  the  line  between  the  aircraft  and  the 
target.  The  calculation  of  this  steering  command  is  the  task  of  the  LOS 
guidance  algorithm.  A  concern  with  the  LOS  guidance  scheme  is  that  some 
missile  control  energy  must  be  used  to  keep  the  missile  on  a  rotating  LOS. 
How  much  control  energy  this  takes  was  not  clear  prior  to  this  effort,  but 
the  missile  control  energy  consumed  is  one  of  the  measures  of 
effectiveness  in  this  study.  This  was  a  concern  because  the  current 
missile  design  turns  with  a  limited  supply  of  pyrotechnic  cartridges.  A 
second  concern  with  this  guidance  method  involved  the  missile's  rocket 
plume  standing  between  the  target  and  the  FLIR  assigned  to  track  the 
target.  This  problem  is  not  addressed  in  this  study,  except  to  note  that 
the  LOS+  guidance  algorithm  was  devised  to  avoid  this  problem.  The 
ability  of  the  FLIR  to  track  a  target  through  the  missile  plume,  cr  to 
detect  a  burned  out  missile  overlaying  a  target  is  a  problem  w.hich 
requires  further  study.  The  LOS  algorithm  analysis  of  Chapter  4  and  the 
conclusions  and  recommendations  of  Chapter  5  address  these  concerns 
further . 

In  working  with  the  LOS  guidance  concept  some  geometry  problems  are 
discovered  that  were  not  discussed  in  any  of  the  literature  uncovered  on 
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Figure  3-1  Rotating  Geometry 


beam  r_ders.  The  algorithm  development  proceeded  as  did  some  of  the  early 
ballistic  missile  algorithms,  with  what  is  called  a  nominal  velocity  and 
tiie  calculation  of  a  velocity  to  be  gained  [7:  65-70], 

The  dot  product  of  this  velocity  vector  to  be  gained  with  the  current 
missile  velocity  vector  provides  the  missile  steering  commands.  The 
nominal  velocity  is  calculated  in  two  parts.  A  missile  that  is  or.  the 
nominal  course,  i.e.  on  the  LOS  between  the  aircraft  and  the  target,  will 
have  a  component  of  velocity  keeping  it  on  the  rotating  LOS,  and  a 
component  of  velocity  directing  it  down  the  LOS  to  the  target.  Tnese 
components  are  shown  in  Figure  3-2  as  and  Vj^.  The  first  component  can 
be  calculated  from  the  LOS  rate  and  the  range  from  the  aircraft  to  the 
missile,  Rm.  The  second  component,  moving  the  missile  down  the  line  of 
sight,  is  simply  the  remaining  missile  velocity  magnitude.  An  additional 
correction  term  must  be  added  to  this  otherwise  nominal  velocity,  to 
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Figure  3-2  Differing  Coordinate  System  of  V_, 


compensate  for  missile  position  and  bring  the  missile  onto  the  target  line 
of  sight.  This  correction  term  is  the  distance  to  the  line  of  sight 
(labeled  position  error  in  Figure  3-1)  divided  by  a  correction  time, 

This  is  the  time  allowed  for  this  position  error  correction  and  it  is  a 
predominate  algorithm  control  parameter.  Thus  the  nominal  velocity  vector 
is  calculated  as  follows: 


V,  t  +  V,.  r 

V,  =  (r,  Ru+  r.  6/tJt  +  (V/-(r,  B, ,)=■)’'"  r  (3-1) 

Where:  =  Nominal  velocity  vector  for  on  trade  missile 

Vj„  =  the  V„  tangent  to  the  line  of  sight 
Vi^  =  the  V„  radial  along  the  line  of  sight 
6  =  B,, 

t^  =  an  input  control  parameter  of  correction  time 
r^  =  horizontal  range  from  aircraft  to  missile 
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B;.  =  LOS  angle  rate  from  aircraft  to  target 
t  =  unit  vector  component  tangent  to  Bjt 

=  Velocity  of  missile  with  respect  to  aircraft 
r  =  unit  vector  component  radial  along 
Note  that  the  correction  for  the  position  error  is  assumed  to  be 
entirely  in  the  tangent  component  of  the  nominal  velocity.  This  is  true 
if  the  difference  angle,  5,  is  small.  Should  this  angle  become  to  large, 
as  in  the  LOS+  guidance  algorithm  discussed  later,  the  portion  of  this 
position  correction  that  lies  in  the  radial  direction  must  also  be 
calculated.  If  the  difference  angle,  5,  grows  large,  a  radial  component 
of  the  correction  term  is  required  to  maintain  this  algorithm  accurately. 
(It  can  be  seen  that  when  6  approaches  90  degrees  a  tangent  term  will  pull 
the  missile  back  towards  the  aircraft  rather  than  towards  the  target. 
This  problem  is  encountered  during  the  LOS+  development  because  a  forecast 
position  is  used  which  caused  8  to  exceed  90  degrees.)  For  a  normal  LOS 
guidance  application,  where  the  target  and  missile  remain  in  the  same 
line-of-sight  angle  quadrant,  this  tangent  component  is  all  the  correction 
term  that  is  necessary.  With  a  nominal  velocity  calculated,  one  must  now 
calculate  the  velocity  of  the  missile  and  translate  it  into  the  same 
reference  frame.  The  missile  velocity  is  calculated  from  the  measurement 
model  parameters  as  follows: 

V,  =  r,  fi,,  t'  +  r„  r'  (3-2) 

Where:  V.  =  missile  velocity  with  respect  to  the  aircraft 
r.  =  horizontal  range  form  aircraft  to  missile 
B,.  =  LOS  angle  rate  from  aircraft  to  missile 
t'  =  unit  vector  component  tangent  to  Bi„ 
r'  =  unit  vector  component  radial  along  B;^ 
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It  can  be  seen  in  Figure  3-2  that  a  rotation  of  the  coordinate  system  is 
required  before  mathematical  operations  can  be  accomplished  on  the  missile 
LOS  vector.  This  vector  rotation  is  accomplished  as  follows: 

=  (r„  cos  (5)  +  r.  sin  (6))  t 

+  (r.  sin  (6)  +  r.  cos  (8)  )  r  (3-3) 

Where:  5  =  the  angle  between  r  and  r' 

In  order  to  understand  the  control  problem  and  construct  the  control 
block  diagram  a  small  angle  approximation  can  be  made  if  the  missile  is 
close  to  its  nominal  position  on  the  aircraft  to  missile  line-of-sight . 
This  approximation  is  accurate  for  angles  of  less  than  10  degrees.  V^ith 
this  assumption  the  missile  velocity  vector  of  equation  (3-3)  reduces  to: 

V.  =  r^  Bj.  t  +  r„  r  (3-4) 

Now  that  the  velocity  of  the  missile  and  the  desired,  or  nominal 
velocity  are  in  the  same  coordinate  system  the  velocity  to  be  gained  (V,) 
is  simply  the  difference: 

V  =  V  -  V 

g  r;  r 

Combining  equations  (3-1)  and  (3-4): 

V,  =  r^  (8/t^+  iJ,,-  Bi„)  t  +  (r,_  -  r„)  r  (3-5) 

Where:  r^  =  (V^^  -  (r„  Bj,)^)’''^ 

Once  the  velocity  to  be  gained  vector  is  computed  the  heading  change 
may  be  computed  by  ta)cing  the  cross  product  as  follows: 

01,.,  =  K3/ivJ  (V,  X  VJ  (3-6) 

0,.,.,  =  K3  r./|vj  (r.  (B.,+8/t,)-  r,  B,J  (3-7) 

This  is  the  equation  used  in  the  LOS  guidance  algorithm.  Notice  that 
the  missile  velocity  magnitude  is  normalized  out  of  the  heading  command 
but  the  velocity  to  be  gained  magnitude  is  not.  The  algorithm  gain 
control,  K3,  must  be  quite  small  to  account  for  this  remaining  magnitude. 
Some  dot  product  steering  methods  normalize  to  one  and  use  K3  to  control 
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all  the  loop  gain.  The  magnitude  of  Vg  was  left  in  the  command  because  it 
represents  a  magnitude  of  required  error  correction. 

Notice  that  these  equations  handle  the  azimuth  guidance  in  the 
horizontal  plane.  An  identical  set  of  these  equations  has  tc  be 
implemented  in  the  vertical  plane  to  calculate  the  missile  pitch  rate 
command,  ©jcmd  • 

Notice  also  that  care  is  taken  to  implement  this  algorithm  using  only 
the  variable  states  available  from  the  measurement  model  described  in 
Paragraph  2.6.  Thus  the  LOS  guidance  algorithm  is  tested  with  identical 
scenarios  available  to  the  other  guidance  algorithms. 

3.2.1  LOS  Algorithm  Control  Parameters.  The  algorithm  allows  two  variable 
inputs,  K3  and  t^,  that  change  the  dynamics  of  the  guidance  algorithm.  By 
studying  the  block  diagram  of  Figure  3-3  it  can  be  seen  that  these 
parameters  control  the  system  gain  and  a  control  loop  zero.  To  derive  the 
block  diagram  and  the  root  locus  plots  requires  some  assumptions  shown  in 
the  following  development. 

Equation  (3-7)  can  be  simplified  by  assuming  that  the  LOS  angle  to  the 
target  and  the  LOS  angle  to  the  missile  will  always  be  close  to  each 
other,  say  within  0.17  radians.  This  is  true  when  the  missile  is  being 
properly  guided  down  the  line-of-sight .  It  is  not  true  during  the  first 
moments  after  launch  for  an  off  angle  shot  (i.e.  an  off  angle  of  greater 
than  0.17  radians)  but  the  missile  quickly  comes  into  this  limitation  so 
the  general  control  principle  being  developed  here  is  effective.  This 
small  angle  assumption  used  to  analyze  the  control  problem  yields  sin(5)= 
6  and  cos(6)=  1.  Assuming  the  missile  is  heading  almost  directly  away 
from  the  aircraft,  a  typical  guidance  situation,  r^  rate  =  r^  rate  =  V. . 
Using  tiiese  approximations  to  simplify  equation  (3-7)  we  see  that: 

r^  B;,  -  sin(©;^  -B,J 

solving  for  IJ,^  and  applying  the  small  angle  assumption  yields: 

B,.  =  V„/r„  (0„-B.,) 


3-7 


Figure  3-3  LOS  Control  Block  Diagram 


VJith  the  velocity  assumption  this  simplifies  equation  (3-7)  into: 

=  K3  r„  (iJu-  Bim+  Bi./tc-  Bu/tc  ) 

With  these  simplifications  in  place  the  block  diagram  of  Figure  3-3  is 
constructed.  The  poles  and  zeros  of  this  control  system  are  plotted,  and 
it  is  seen  that  t^.  controls  a  z  'ro  placement.  The  plotting  of  the  root 
locus  here  is  not  a  trivial  matter  because  there  is  a  pole  location  at 
-V^/r^.  As  the  missile  accelerates  and  leaves  the  aircraft  this  pole 
location  moves  and  causes  the  root  locus  plot  to  have  a  time  dependence. 
V,/r„  changes  according  to  Table  3-1.  As  the  missile  gets  farther  away 
from  the  aircraft  the  LOS  guidance  control  characteristics  change.  Three 
snapshots  of  the  root  locus  plots  are  shown  in  Figure  3-4.  The  three  plots 
represent  times  of  .06,  .5,  and  1.0  seconds,  with  a  zero  set  by  t„  at  -10. 
These  plots  are  representative  of  the  root  locus  variations  during  a 
launch  with  a  LOS  guidance  scheme.  They  reveal  the  time  variation  of  this 
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Table  3-1  Typical  flight  V-/R-  pole  values 


Time 

sec 

V. 

ft /sec 

K 

k  ft 

-  V„/R„, 
/sec 

0.02 

55 

0.005 

-11.83 

0.04 

110 

0.006 

-17.80 

0.06 

166 

0.009 

-18 . 69 

0.08 

227 

0.013 

-17.28 

0 . 10 

284 

0.018 

-15.52 

0 . 12 

340 

0.025 

-13 . 87 

0 .14 

398 

0.032 

-12.44 

0 .16 

455 

0.041 

-11.23 

0.18 

508 

0.049 

-10.30 

0.20 

567 

0.060 

-9.43 

0.50 

1562 

0.368 

-4.24 

1.00 

3322 

1.574 

-2.11 

1.50 

5565 

3.777 

-1.47 

2 . 00 

5123 

6.445 

-0.795 

2.50 

4720 

8.892 

-0.531 

2.89 

4431 

10.66 

-0.41 

root  locus  throughout  the  flight.  A  more  efficient  control  could  be 
developed  by  tuning  the  system  gain  for  various  phases  of  the  flight.  The 
gain  could  then  be  adjusted  to  provide  the  desired  damping  for  the  entire 
flight.  In  this  development  a  constant  gain  was  used  to  provide  favorable 
damping  over  the  range  of  control  characteristics  found  in  a  typical 
flight  of  the  defined  missile.  The  placement  of  the  zero  and  the  gain  is 
controlled  by  varying  t^  and  K3 .  This  is  accomplished  to  give  best 
algorithm  performance  and  is  the  first  thing  explored  in  the  analysis 
section  of  Chapter  4. 

3 . 3  Line-of-siqht  Plus  (LOS-*-)  Algorithm  Development 

The  LOS+  algorithm  is  a  unique  solution  devised  to  move  the  missile  off 
the  target  line-of-sight  (in  an  attempt  to  ease  the  FLIR  target  and 
missile  discrimination)  and  to  reduce  the  control  energy  involved  in 
following  a  rotating  line-of-sight  [1].  The  LOS+  guidance  algorithm  is 
not  new  but  no  references  were  found  that  give  its  description  or  document 
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Figure  3-4  Root  Locus  Plots  During  a  Typical  Flight 


3-10 


its  use.  As  with  the  LOS  guidance,  however,  the  application  of  the 
algorithm,  with  a  moving  tracker,  is  unique.  The  principle  involved  in 
this  guidance  concept  requires  forecasting  the  aircraft  position  and  the 
LOS  to  the  target,  Bi-  and  ^t  the  time  of  missile  intercept  and  flying 
the  missile  down  this  reasonably  stable  line-of-sight  into  the  target . 
Figure  3-5  is  a  layout  of  this  concept. 

Two  challenges  evident  in  this  concept  are  the  requirement  to  calculate 
the  time  to  missile  impact  and  the  control  problem  involved  in  guiding  the 
missile  on  to  this  future  line-of-sight.  In  the  very  early  development  of 
this  guidance  algorithm,  its  inefficiencies  began  to  show  through.  The 
first  evident  drawback  is  the  missiles  high  velocity  when  the  first  major 
turn  onto  the  future  line-of-sight  is  encountered.  This  high  velocity 
causes  this  turn  to  impose  a  high  control  energy  requirement  as  seen  in 
the  analysis.  Figure  3-5  shows  this  problem,  though  somewhat  exaggerated 
because  of  the  high  off  angle  shot  pictured. 

Attempts  were  made  to  "aim"  for  different  intercept  points  on  the 
future  line-of-sight  but  each  solution  is  a  compromise  in  either  accuracy 
or  control  energy.  Because  of  its  high  velocity  and  short  flight  time  an 
optimal  solution  must  be  developed  to  efficiently  put  the  missile  on  the 
future  line-of-sight.  Since  the  optimal  solution  will  involve  a  time 
varying  gain  and  Locus  pole/zero  placement  (reference  the  LOS  control 
parameters  of  Para  3.2.1)  it  is  considered  beyond  the  scope  of  this 
research.  The  other  guidance  algorithms  developed  here  were  considered 
more  advantageous  than  spending  much  effort  on  this  algorithm  development. 

A  second  drawback  to  the  LOS+  algorithm  is  its  sensitivity  to  the 
aircraft's  maneuvers  after  launch.  Yogi  Bera  once  said  "the  future  is  not 
what  it  used  to  be,"  and  this  is  especially  true  of  a  future  line-of- 
sight  when  a  pilot  pulls  an  unpredicted  5  G  turn.  The  other  algorithms 
developed  here  are  intuitively  more  independent  of  the  launch  aircraft's 
maneuvers.  Again  the  LOS+  algorithm  could  be  better  tuned  to  account  for 
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aircraft  acceleration  but  the  work  involved  is  deemed  more  tedious  than 
advantageous,  given  the  shortfalls  of  this  algorithm.  The  challenges  of 
calculating  an  impact  time  and  guiding  to  the  future  line-of-sight  are 
discussed  in  the  following  paragraphs. 

3.3.1  Missile  Impact  Time  Calculation.  The  missile  impact  time  must  be 
calculated  in  order  to  project  the  aircraft  position  properly.  The  missile 
impact  time  is  calculated  with  a  simple  binary  search  iteration  technique. 
The  time  to  go,  t^  is  calculated  as  follows: 

t,  =  P.,  /  |V„,|  (3-8) 
Where  R,  =  range  to  go  i.e.  missile  to  target  range  (see  equation  (3-9)) 

I  V„.  I  =  Average  missile  Velocity  (see  equation  (3-10)) 

A  simple  conversion  to  the  NED  x,  y,  and  z  coordinates  and  a  magnitude 
calculation  of  their  differences  is  used  in  this  calculation.  This  is 
accomplished  using  the  following  equations: 
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R,  =  [(r^  sin  (Bit)  -  sin(6i„))'+ 

(rt  cos  (Bit)  “  cos(Bi„))^+ 

(Rt  sin(B2t)  -  Rr,  sin(Bj„)  )^]’ 


(3-9) 


The  average  missile  velocity  is  calculated  by  a  weighted  time  average 
of  the  acting  accelerations  as  follows: 


IF  t  <  tp 


lV„t|  =  V„  +  1/2  (at,  (t. 


t  )  +  a„  t.  ] 


IF  t  >=  t. 


V„t  =  V„  +  1/2  a,t  t. 


(3-10) 


(3-11) 


Where:  Vj  =  missile  velocity  at  time  of  calculation  (see  equation  (3-14)) 
t  =  time  of  calculation 
tg  =  time  to  go  before  impact 
teuRN  “  time  of  rocket  burn 

at,  =  Average  accel  caused  by  missile  thrust  (see  equation  (3-15) ) 
a<,t  =  Average  accel  caused  by  missile  drag  (see  equation  (3-16)) 

The  present  missile  velocity  is:  =  |v^  -  Vt( 

To  perform  this  subtraction  the  target  velocity  term  is  rotated  into 
the  line-of-sight  to  the  missile  as  follows: 


^Im  ■*"  Rri  ^3  (3-12) 

Vt  =  [r,  cos  (5)  +  rt  B;t  sin(6)]  1,  + 

[rt  sin  (5)  +  r,  Bi,  cos  (6)]  I,  + 

Rt  Bjt  I3  (3-13) 

Subtracting  the  terms  and  talcing  the  magnitude  yields: 

V(,  =  (  (  r,-  (rt  cos  (6)  +  r,  B,.  sin (6))  )'  + 

(  r„  Bi„  -  (-r.  sin  (6)  +  rt  B,t  cos  (5))  )^  + 

(R^  B^t,-  Rt  fljt)^  (3-14) 

Where:  r„ =  horizontal  range  from  aircraft  to  missile 

r,  =  horizontal  range  rate  from  aircraft  to  missile 
B,„  =  LOS  angle  from  aircraft  to  missile 


B,„  =  LOS  angle  rate  from  aircraft  to  missile 


r.  =  horizontal  range  from  aircraft  to  target 
r.  =  horizontal  range  rate  from  aircraft  to  target 
=  LOS  angle  from  aircraft  to  target 

iii.  =  LOS  angle  rate  from  aircraft  to  target 
=  range  form  aircraft  to  missile 

jij„  =  Pitch  LOS  angle  rate  from  aircraft  to  missile 
Rt  =  range  form  aircraft  to  target 

iij.  =  Pitch  LOS  angle  rate  from  aircraft  to  target 

The  acceleration  and  drag  terms  required  for  equation  (3-10)  and  (3-11) 
are  averages  calculated  from  missile  input  parameters.  The  use  of  these 
average  accelerations  over  the  remaining  time  of  flight  meets  the  accuracy 
requirements  of  this  simple  time  to  go  calculation.  The  acceleration  and 
drag  terms  are  calculated  as  follows: 

a,^  =  (TH)  *  G  /  (Wtbo  +  1/2Wtpu)  (3-15) 

AND: 

-a  JC  D'  V.  (V,  +  l/2a,,  teuBN) 

Sdr  =  - 

8  *  WTBD/G  (3-16) 

Where:  TH  =  Rocket  Thrust  input  from  missile  data  file 
G  =  Gravitational  constant 

Wtbd  “  missile  body  weight  input  from  missile  data  file 
Wtfu  =  missile  fuel  weight  input  from  missile  data  files 
o  =  air  density  (.0023  slugs/ft’) 

D  =  missile  diameter  input  from  missile  data  files 
Vj  =  Speed  of  Sound  at  average  altitude,  1087  ft/sec 
Vj  =  aircraft  velocity 

These  equations  are  now  used  in  a  binary  search  iteration  to  calculate 
the  time  to  go,  t,.  The  binary  search  is  implemented  as  follows: 

Step  1 :  Calc  R,  and  V^,  and  guess  at  time  to  go, 

Step  2:  Calculate  |v„J  via  equation  (3-10)  or  (3-11) 

Step  3:  Calculate  t,  via  equation  (3-8) 

Step  4:  Set  t,,.„  =  t,„„  +  1/2  (t,+  t,„„,) 
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step  5:  IF  ABS  >  Tolerance  Repeat  from  Step  2 

This  method  converged  on  to  the  missile  time  to  go  value  in  about  10 
iterations  in  each  case  tested.  It  also  makes  the  assumption  that  the 
missile  is  using  all  its  velocity  to  get  straight  to  the  target. 

3.3.2  LOS-t-  Control  Challenge.  Getting  a  missile  of  this  speed  to  fly 
down  a  predicted  aircraf t-to-target  LOS  in  the  short  time  allowed,  turns 
out  to  be  a  challenge.  The  approach  taken  here  is  one  of  several  that 
could  work,  but  others  were  not  attempted  because  there  are  more  efficient 
optimal  control  solutions  that  should  be  used  if  this  guidance  methodology 
were  to  be  considered  further. 

This  approach  follows  the  footsteps  of  the  previous  LOS  development 
that  used  a  nominal  velocity  calculation,  a  velocity  to  be  gained,  and  a 
dot  product  steering  command  calculation.  The  complete  nominal  velocity 
calculation  involves  two  factors.  A  nominal  velocity  for  an  on  track 
missile  and  a  position  correction  term  are  required  to  move  the  missile 
towards  this  'on  track'  condition.  The  nominal  velocity  for  an  on  track 
missile  is  simply  a  velocity  moving  the  missile  down  the  stable  future  LOS 
to  the  target.  The  nominal  velocity  magnitude  is  set  equal  to  the  current 
missile  velocity  magnitude.  So  its  calculation  requires  only  a  rotation 
of  the  missile  velocity  onto  the  future  LOS.  This  rotation  is 
accomplished  as  follows.  The  future  LOS  is  calculated  assuming  a  constant 
2;.  rate: 

li',.  =  B...  +  B,,  t. 

The  horizontal  missile  velocity  magnitude  is  calculated: 
v„  =  [r/  +  (r,  BiJ'"  ]"'■ 

and  the  velocity  is  rotated  into  this  velocity  frame: 

V..,  =  V,  cos  (fli.-B' ,, )  r  +  v„  sin  (B.„-B' ,, )  t  (3-17) 

where:  B'  .  =  Future  LOS  to  target  (from  North) 

V.  =  Horizontal  missile  velocity  magnitude 

=  Nominal  velocity  without  position  correction 
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r  =  unit  vector  along  missile  velocity  path 
t  =  unit  vector  tangent  to  missile  velocity  path 
There  are  several  ways  to  calculate  the  position  error.  They  depend  on 
where  the  missile  is  projected  on  the  future  LOS.  The  selected  scenario 
based  on  performance,  is  shown  in  Figure  3-6.  The  correction  distance,  1^, 
is  calculated  as  follows: 

R,  =  V.,  t,  -d 

where:  d  =  R„  cos  (©i,  -  Ili„)  -  a 

and:  tan(©i,  -  li'it)  =  R„  sinO^  -  film) /a 

Combining  these  terms  yields: 

Rc  =  v,ct,-R„[cos  (©la-fiim)  -  sin(0i,-fii„) /tan(©i,-fi' u)  ] 

As  shown  in  the  figure  the  direction  of  Re  is  chosen  parallel  to  the 
aircraft  heading.  Expressing  R^  as  a  vector  in  the  fii„  frame  and  dividing 
by  the  time  allowed  for  the  position  correction  yields: 

=  Re/t<.  [cos(0,.-  film)  r  +  sin{©,^-  fiim)t]  (3-18) 

where:  V^,  =  position  correction  portion  of  nominal  velocity 
The  total  nominal  velocity  is  now: 

V  =  V  +  V 

''n  ''na  ''nb 

Vn  =  |Vm|  cos  (film  -  fi'u)  +  Rc/tc  COS(©.^  -  fijm)  t 

+  |V„|  sin(Bi„  -  fi'iJ  +  Re/te  sinOi,  -  &!„)  t  (3-19) 

The  velocity  to  be  gained  is  and  the  commanded  missile  heading 

rate  is  proportional  to  the  dot  product  of  Vm  with  as  follows: 

0^,  =  K3  (V,  *  VJ  (3-20) 

To  ma)te  this  algorithm  perform  better  it  is  necessary  to  weight  the 
position  correction  term  more  heavily  during  the  final  stage  of  the 
engagement  than  at  the  beginning.  This  is  not  attempted  in  this  research 
because  this  algorithm  is  inferior  to  the  LOS  and  Proportional  guidance 
algorithms . 
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Notice  that  these  equations  handle  the  azimuth  guidance  in  the 
horizontal  plane.  An  identical  set  of  equations  must  be  implemented  in 
the  vertical  plane  to  calculate  the  missile  pitch  rate  command, 

Care  is  taken  to  implement  this  algorithm  using  only  the  variable 
states  available  from  the  measurement  model  described  in  Paragraph  2.6. 
Thus  the  LOS+  guidance  algorithm  is  tested  with  identical  scenarios 
available  to  the  other  guidance  algorithms. 

3 . 4  Pursuit  Guidance  Algorithm  Development 

The  pursuit  guidance  algorithm  is  an  old  and  simple  guidance  solution 
whereby  the  missile  is  continually  flown  at  the  known  target  position. 
This  is  a  relatively  easy  implementation  to  use  and  the  geometry  is  shown 
in  Figure  3-7.  The  algorithm  requires  the  calculation  of  the  missile  to 
target  line-of-sight  and  the  missile  heading,  or  at  least  its  velocity 
vector.  These  two  calculations  are  accomplished  as  follows: 
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From  Figure  3-7  it  is  seen  that: 

r  r^  sin(8) 


a. 


tan"' 


r^  cos  (8)  -  r„ 


1 

J 


=  tan-‘(r,  iii„/  r„) 
where:  8  =  C..- 

Then  the  commanded  azimuth  turn  is  made  proportional  to  a- 
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(3-22) 


+  ttj  as  follows : 


© 


ICfT.d 


-Kj  (a,  +  a,  ) 


(3-23) 


Notice  that  these  equations  handle  the  azimuth  guidance  in  the 
horizontal  plane.  An  identical  set  of  equations  are  implemented  in  the 
vertical  plane  to  calculate  the  missile  pitch  rate  command,  0,cna  ■  The 
Pursuit  guidance  algorithm  is  a  straight-forward  implementation.  Again 
notice  that  care  is  talcen  to  implement  this  algorithm  using  only  variables 
available  from  the  measurement  model  described  in  paragraph  2.6.  Thus  the 
pursuit  guidance  algorithm  is  tested  with  the  same  scenarios  and  data 
calculations  as  the  other  algorithms  under  test. 


3-18 


J  .  D 


Proportional  Guidance  Algorithm  Development 


The  proportional  guidance  algorithm  uses  the  line-of-sight  rate  as  a 
guidance  parameter  instead  of  the  line-of-sight  angle  used  in  the  pursuit 
guidance  [8:198]  .  The  geometry  of  the  problem  is  identical  to  the  pursuit 
geometry  layout  shown  in  Figure  3-7  above.  Here  the  guidance  algorithm 
commands  a  missile  heading  rate-of-change  which  is  proportional  to  the 
rate-of-change  of  the  line-of-sight  from  the  missile  to  the  target .  The 
proportional  constant/  K3,  is  called  the  navigation  constant  [8:199].  A 
computational  challenge  is  the  need  to  calculate  the  missile-to-targct 
line-of-sight  rate  .  In  the  normal  implementation  of  this  guidance 
algorithm  Cij„.  rate  is  simply  measured  from  the  missile's  target  tracker. 
In  this  application  this  value  is  constructed  from  the  measured  missile 
velocity  and  target  velocity.  Thus  the  proportional  guidance  algorithm  is 
implemented  as  follows: 


©icmd  ~  (3  ^4) 

where:  0icn,j  =  the  missile  azimuth  turn  command 
Kj  =  the  navigation  constant 

=  the  missile  to  target  line-of-sight  rate 
rate  remains  to  be  calculated.  The  calculation  of  the  target  velocity 
relative  to  the  missile  is  accomplished  using  the  V*  and  Vt  as  expressed 
in  equation  (3-12)  and  (3-13)  .  It  follows  then,  that  the  target  velocity 
relative  to  the  missile  is: 


=  V,  -  V,  =  V,  1.  +  V;  1, 


(3-25) 


Where:  V.  =  r„  -  r.  cosdi...  -  B,„)  -  r^  Ii,t  sindi,..  -  lii,) 

V;  =  r^  li;.  -  r.  sin(li,.  -  Bi„)  -  r.  ii..  cosdi:-  -  Iii„.) 

The  component  of  this  relative  velocity  that  causes  a  rotation  of  the 
missile-to-target  line-of-sight  is  tangent  to  the  missile-to-target  line- 
of-sight,  liip.  .  This  tangent  velocity  term  can  be  expressed  two  ways,  as 
follows : 
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film'.  =  ^2  cos(e.i^.  -  +  Vj  sin(Ii,^  -  13-26) 

Solving  these  two  for  yields: 

Vj  cos(Bi„-  Bi„,)+  V,  sin(B,„-  Bi„J 

B:,.  =  - 

r..t  (3-27) 

The  Bi„,.  angle  used  in  this  equation  is  calculated  from  the  NED  X  and  Y 
coordinates  as  follows: 


r,  =  horizontal  range,  aircraft  to  missile 
r^  =  horizontal  range  rate,  aircraft  to  missile 
B;„  =  azimuth  line-of-sight,  aircraft  to  missile 
B.,  =  azimuth  line-of-sight  rate,  aircraft  to  missile 
Notice  that  these  equations  handle  the  azimuth  guidance  in  the  horizontal 
plane.  An  identical  set  of  these  equations  are  implemented  i.'  the 
vertical  plane  to  calculate  the  missile  pitch  rate  command,  • 

Again  care  is  ta)cen  to  implement  this  algorithm  using  only  the  variable 
states  available  from  the  measurement  model  described  in  paragraph  2.6. 
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Thus  the  proportional  guidance  algorithm  is  tested  with  the  same 
scenarios  and  data  calculations  as  the  other  algorithms  under  test . 


3 . 6  Range  and  Range  Rate  Estimator 

In  the  event  that  the  aircraft  will  launch  without  electromagnetic 
emittions  it  is  necessary  to  test  each  algorithm's  performance  with  a 
range  estimator  rather  than  range  and  range  rate  from  a  radar.  To 
accomplish  this  test,  a  very  simple  range  and  velocity  estimator  is  built 
to  use  in  place  of  the  radar  measurements.  To  do  this  average  thrust  and 
drag,  calculated  in  equation  (3-15)  and  (3-16),  is  used.  An  Euler 
integration  is  used  to  calculate  the  velocity  and  another  to  calculate  the 
range.  This  rough  estimator  assumes  the  missile  is  always  accelerating 
directly  away  from  the  launch  aircraft,  an  assumption  that  is  accurate 
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Figure  3-8  Range  Estimator  Performance  Comparison 

except  for  high  off-angle  shots.  A  graph  of  this  range  estimation  is 


shown  in  Figure  3-8.  The  actual  range  shown  in  this  graph  is  taken  from  a 


LOS  guidance  shot  with  a  30  degree  off  angle  initial  heading.  The  range 
rate  estimation  is  made  by  taking  the  estimated  missile  velocity  and 
backing  out  the  components  causing  the  azimuth  rate  and  pitch  rate.  This 
is  done  because  a  missile  azimuth  and  pitch  rate  can  be  obtained  from  the 
FLIR.  Thus  the  range  rate  is  calculated  as  follows: 

R.  =  R„t  cos(Ii,J]^  -  (3-30) 


This  range-rate  estimate  as  well  as  the  actual  range-rate  from  the  LOS 
guidance  shot  is  shown  in  Figure  3-9.  Although  more  accurate  on-board 
estimates  of  missile  range  and  range  rate  could  easily  be  developed,  this 
one  is  used  to  demonstrate  some  of  the  algorithms  dependence  on  range  and 
range  rate.  The  objective  here  is  not  to  build  an  accurate  estimator, 
just  a  consistent  one  that  can  be  used  in  each  test. 
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IV .  Algorithm  Analysis 


4 . 1  Introduction 

The  performance  of  each  algorithm  developed  in  chapter  3  is  tested  and 
the  results  are  explained  in  this  section.  In  performing  analysis  there 
are  always  more  scenarios  desired  than  one  can  possibly  evaluate.  For  the 
initial  analysis  a  30  degree  off-angle  shot  is  chosen  because  it  strains 
the  algorithms  initial  correction  capabilities.  This  is  used  to  test  and 
establish  the  gain  and/or  control  parameters  for  each  algorithm.  For 
this  shot  uhe  target  is  at  030  degrees  LOS  with  a  heading  of  due  east, 
(i.e.  090  degrees).  The  aircraft  and  missile  initial  heading  is  060 

degrees.  Then  a  simpler  5  degree  off-angle  scenario  is  used  to  test 
reactions  to  an  aircraft  5-G  turn  during  the  missile  flight.  Here  the 
target  is  at  055  degrees  LOS  with  a  heading  of  due  west,  (i.e.  270 

degrees).  The  aircraft  and  missile  initial  heading  is  060  degrees. 
Finally  the  aircraft's  radar  is  turned  off  and  a  simple  missile  range  and 
velocity  estimator  is  used  to  'sense'  the  missile  range  and  range  rate. 
This  process  shows  the  error  in  the  algorithms  when  accurate  range  and 
range  rate  are  not  available  to  the  algorithm.  The  "radar  off"  test 
repeated  the  30  degree  off-angle  scenario.  The  only  algorithm  that  has 
performance  reasonably  independent  of  these  range  parameters  is  the  LOS 
algorithm.  This  chapter  is  organized  to  analyze  the  data  for  each  of 
these  algorithm  tests . 

4 . 2  Figures  of  Merit 

The  missile  miss  distance  and  the  amount  of  expended  control  force  is 
used  as  the  primary  figures  of  merit  (FOM)  for  this  analysis.  Each 
algorithm  can  be  adjusted  from  an  under-damped  to  an  over-damped  condition 
with  the  control  parameters.  In  general  this  changed  the  miss  distance 
from  small  to  large  and  the  control  energy  from  large  to  small.  So  by 
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observing  the  primary  FOMs  a  cont  :ol  parameter  can  be  chosen  to  obtain  the 
best  algorithm  performance.  Thus,  minimizing  these  two  FOMs  forms  the 
basis  for  most  of  this  analysis. 

The  missile  miss  distance  is  measured  during  the  missile  flight  by 
recording  the  missiles  closest  approach  to  the  target .  There  is  no 
consideration  given  for  a  missile  impacting  the  ground  before  it  reaches 
this  closest  approach.  One  should  exercise  caution  in  using  this  FOM  to 
obtain  a  probability  of  kill.  The  miss  distance  FOM  is  a  tool  to  evaJ.uate 
the  guidance  algorithm  accuracy. 

The  amount  of  control  force  required  to  steer  the  missile  is  the  FOM 
used  to  measure  the  guidance  algorithms  efficiency.  This  control  force  is 
measured  by  integrating  the  commanded  turn  rates  throughout  the  missile 
flight.  This  is  accomplished  with  the  following  equation: 

E„  =  E„  +  |vj  M  dt 

Where:  £„  =  total  control  moment  commanded  (  lbs  sec) 

V„  =  total  missile  velocity  (ft/sec) 

©land  “  commanded  missile  heading  rate  (rad/sec) 

©2aBd  =  commanded  missile  pitch  rate  (rad/sec) 

M  =  missile  mass  (slugs) 

dt  =  integration  time  of  one  execution  cycle  (sec) 

Although  this  symbology  hints  that  this  is  a  measure  of  the  control 
energy,  it  is  not.  Monitored  here  is  the  time  integral  of  turning  force 
commanded,  integrated  over  the  time  of  the  flight.  To  calculate  the 
control  energy  requires  more  information  on  how  efficiently  the  missile 
transforms  the  commands  to  turning  motion,  which  requires  much  more  detail 
than  is  available  about  the  missile.  The  measure  of  the  total  force  and 
time  product  that  it  is  commanded  gives  a  force  time  unit  (FTU)  measure 
that  is  useful  in  comparing  the  algorithms  even  if  it  is  not  physically 
realizable.  The  control  FTU  is  a  useful  FOM  for  determining  the  guidance 
algorithms  efficiency. 


4-2 


4 . 3  LOS  Control  Parameter  Testing 


The  various  runs  used  to  establish  the  LOS  algorithm  control  parameters 
K3  and  t(,/  are  shown  in  Figure  4-1  through  4-5.  The  only  parameters 
changed  for  these  runs  are  the  control  parameters .  Included  in  these 
figures  is  a  scenario  ground  map  that  pictures  the  aircraft  track,  the 
target  track,  and  the  missile  flight  path.  In  the  upper  left  corner  of 
this  map  is  the  type  of  algorithm  used,  some  figures  of  merit  and  the 
control  parameters.  Inset  into  each  of  the  scenario  maps  is  a  graph  of 
the  heading  rate  commands  and  a  graph  of  the  missile  heading  and  pitch. 

The  graph  of  the  missile  heading  shows  the  effect  of  the  LOS  control 
variables  t^  and  Kj.  In  Figure  4-1  notice  the  slight  heading  overshoot 
around  the  0.8  second  TOF  and  the  quick  settling  on  to  a  steady  heading 
rate  that  effectively  brings  the  missile  into  the  target.  Notice  also  the 
change  of  heading  from  about  32  degrees  to  about  28  degrees  during  the 
last  two  seconds  of  flight.  This  is  caused  by  the  rotating  aircraft  to 
target  line  of  sight  as  explained  in  chapter  three. 

A  change  in  the  heading  command  is  evident  at  the  1.5  sec  TOF,  the  time 
that  the  missile  rocket  engine  burns  out.  This  is  caused  by  the  step 
change  in  the  missile  acceleration  that  occurs  at  this  point.  Recall  that 
the  missile  accelerations  were  not  accounted  for  in  the  algorithm 
development,  but  it  does  effect  the  missile  LOS  and  LOS  rate.  This  causes 
the  control  system  to  have  an  increase  in  accuracy  when  the  unaccounted 
for  acceleration  just  decreased  drastically.  This  step  through  the 
control  system  settles  quickly  and  is  barely  discernable  in  the  heading 
plot . 

Figure  4-2  through  4-5  show  various  stages  of  control  damping  from 
under-damped  to  over-damped.  Recall  from  chapter  three,  that  a  control 
zero  is  at  -  l/t^  and  the  gain  varies  with  Kj,  and  that  a  pole,  located 
at  -V„  /R„,  moves  during  the  missile  flight.  A  typical  pole  location  for 
several  TOF  values  is  shown  in  Table  3-1.  This  "movement"  of  the  control 
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Figure  4-2  LOS  Control  Parameter  Test  #2 
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Figure  4-3  LOS  Control  Parameter  Test  #3 
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Figure  4-4  LOS  Control  Parameter  Test  #4 
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Figure  4-5  LOS  Control  Parameter  Test  #5 


pole  location  causes  the  control  characteristics  to  change  during  the 
missile  flight.  It  is  .very  difficult  to  see  this  change  in  control 
characteristic  in  the  figures,  but  the  objective  is  an  average  control 
response  that  acts  critically  damped  throughout  the  flight. 

Figure  4-6  shows  all  the  FOM  results  for  the  LOS  algorithm  runs.  The 
first  five  graphs  show  the  tests  run  for  control  parameter  selection.  By 
analyzing  the  results  it  is  apparent  that  the  best  shot  is  obtained  with 
t,  =  0.1  and  K3  =  0.004.  However,  the  best  overall  performer  is  t^  =  0.1 
and  Kj  =  0.003.  Notice  from  the  control  graphs  that  these  control 
parameters  give  smoother  damped  sinusoidal  control  curves.  In  running 
with  Kj  =  0.004  the  control  is  more  erratic  and  did  not  perform  well  in  all 
the  other  scenarios .  (Additional  scenario  runs  are  presented  in  Appendix 
A.)  The  runs  with  Kj  =  0.003  show  the  good  FOM  values  as  well  as  the  best 
damped  control  responses.  The  remainder  of  the  LOS  guidance  tests  are  run 
using  =  0.1  and  Kj  =  0.003. 

4 . 4  LOS  5  Degree  Off-Angle  Tests 

The  less  dynamic  5  degree  off-angle  shot  is  run  to  test  the  algorithm 
with  less  stringent  initial  control  requirements,  and  to  test  the  effect 
of  a  5-G  aircraft  turn  on  the  guidance  situation.  The  bar  graph  labeled 
Sen  #1  and  Sen  #2  in  Figure  4-6  represent  these  two  tests.  Figure  4-7  and 
4-8  give  the  ground  map  and  control  graphs  of  these  two  runs.  In  the  non¬ 
maneuvering  shot  of  Figure  4-7  a  slight  overshoot  and  very  quick  settling 
onto  an  almost  constant  heading  to  the  target  are  evident. 

Figure  4-8  illustrates  the  heading  changes  more  dramatically  as  the 
turning  aircraft  causes  a  change  in  the  LOS  from  the  aircraft  to  the 
target.  The  missile  coiranands  are  a  lot  less  stable  in  this  shot,  in  fact 
the  step  from  the  engine  burnout  at  the  1.5  sec  TOF  point  is  more 
noticeably  pronounced.  The  algorithm  is  strained  here  to  accoun"  for  a 
change  in  the  aircraft  to  target  LOS  rate  that  is  not  accounted  for  in  its 
development.  Thus  the  position  correction  portion  of  the  algorithm  is 
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Figure  4-6  FOM  for  All  LOS  Guidance  Runs 

exercised  to  make  corrections  that  are  not  considered  in  the  algorithm 
development,  where  accelerations  are  not  figured  into  the  LOS  rate 
calculations.  Notice  that  the  accuracy  decreases  in  this  shot  for  the  same 
reason . 

Aircraft  acceleration  could  be  accounted  for  in  the  algorithm  since  it 
is  measured  on  board.  Acceleration  is  not  incorporated  into  the  LOS 
algorithm  because  it  is  not  predictable  and  because  its  effect  can  be 
corrected  via  the  position  correction  portion  of  this  algorithm.  The 
requirement  to  further  tune  the  algorithm  to  account  for  aircraft 
acceleration  may  depend  on  the  algorithm's  sensitivity  to  sensor  errors 
which  have  not  yet  been  explored. 


LOS  Tc/IC3*100D 
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The  performance  of  the  LOS  algorithm  with  critical  damping  and  no 
sensor  error  appears  to  be  superb  in  this  benign  engagement .  However  with 
the  introduction  of  an  aircraft  acceleration  the  miss  distance  doubles. 
This  emphasizes  that  a  hyper-velocity  missile  strains  the  very  best  of  the 
guidance  algorithms. 

4 . 5  Radar  Off  LOS  Tests 

When  the  sensor  model  uses  a  very  simple  3DOF  missile  model  to  estimate 
the  range  and  range  rate  measurements,  the  algorithm's  dependence  on 
accurate  range  data  is  examined.  Figure  4-9  shows  the  no-radar  scenario 
and  it  should  be  compared  to  Figure  4-1  which  shows  the  same  shot  with 
accurate  range  and  range-rate  data. 

Although  there  is  an  increase  in  the  miss  distance  the  algorithms 
performance  is  relatively  unaffected  by  poor  range  and  range  rate  data. 
^  jor  range  data  causes  the  LOS  position  calculations  to  have  a 
proportional  error  but  the  rest  of  the  algorithm  is  independent  of  the 
radar  range.  Recall  from  chapter  3  that  the  missile  range  from  the 
aircraft  does  not  affect  algorithm  performance,  the  algorithms  only 
purpose  is  to  keep  it  on  the  LOS  between  the  aircraft  and  the  target . 
Since  obtaining  an  accurate  radar  measurement  on  this  very  fast  and  small 
missile  is  a  concern,  it  is  a  comfort  that  the  LOS  algorithm  is  not 
dependent  on  an  accurate  radar.  In  fact  it  can  operate  effectively 
without  a  radar. 

4  .  6  I-OSt  Control  Parameter  Testing 

The  various  runs  used  to  establish  the  LOS+  algorithm  control 
parameters  K,  and  t^,  are  shown  in  Figure  4-10  through  4-13.  The  only 
parameters  changed  for  these  runs  are  the  control  parameters.  As  before 
a  scenario  ground  map  is  included  that  pictures  the  aircraft  track,  the 
target  track  and  the  missile  flight  path.  In  the  upper  left  corner  of 
this  map  is  the  type  of  algorithm  used,  some  figures  of  merit,  and  the 
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control  parameters.  Inset  into  each  of  the  scenario  maps  is  a  graph  of 
the  heading  rate  commands  and  graph  of  the  missile  heading  and  pitch. 

The  conglomerations  of  flight  paths  and  control  graphs  found  in  these 
four  figures  only  begin  to  illustrate  the  difficulty  in  tuning  and 
developing  this  algorithm.  Figure  4-10  illustrates  that  the  best 
accuracy  did  not  coincide  with  a  very  efficient  flight  path.  In 
progressing  toward  a  better  flight  path,  poor  control  of  the  position 
error  gives  unacceptable  accuracy.  Since  the  algorithm  is  proposed 
partly  to  reduce  wasted  control  energy,  effort  is  given  to  find  the 
control  values  that  minimize  the  control  force  FOM.  This  is  shown  in 
Figure  4-13  but  notice  that  the  control  force  FOM  doubles  those  used  in 
the  LOS  algorithm. 

The  effort  involved  trying  to  guide  this  missile  to  a  future  line  of 
sight  and  balance  its  efficiency  and  accuracy.  The  sensitivities  caused 
dynamic  fine  tuning  of  the  control  parameters  to  be  required  for  this 
approach.  If  it  is  still  deemed  necessary  to  keep  the  missile  off  the 
target  line  of  sight,  a  variable  gain  will  be  required. 

Examining  Figure  4-13,  the  shot  that  flew  most  efficiently  (but  missed 
by  the  most)  one  can  see  the  sharp  turn  onto  the  predicted  LOS.  The 
commanded  control  then  causes  an  overshoot  and  quick  settling  onto  the 
correct  heading  at  about  the  2  second  TOF  point.  At  about  1/2  seconds 
before  impact  one  see  the  commanded  heading  rate  increase.  This  is  where 
the  time-to-go  prediction  is  discontinued  and  the  algorithm  strives  to 
keep  the  missile  on  the  current  target  LOS.  Recall  that  the  time  to  go 
calculation  is  abandoned  when  t,  <  0.5  sec.  This  is  to  keep  the  t^ 
calculation  errors  from  effecting  accuracy.  (In  hindsight  the  author 
feels  that  it  will  not  effect  it  greatly.)  To  improve  accuracy  the 
control  parameters  need  to  be  altered  during  the  final  engagement.  The 
best  method  of  doing  this  might  involve  an  optimal  control  design. 
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Figure  4-10  I.OS+  Control  Parameter  Test  #1 
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Figure  4-11  LOS+  Control  Parameter  Test  #2 
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Figure  4-13  LOS+  Control  Parameter  Test  #4 


Figure  4-14  shows  all  the  FOM  results  for  the  LOS+  runs.  The  first 
four  graphs  show  the  tests  run  for  control  parameter  selection.  By 
analyzing  the  results  of  the  first  four  runs  it  is  apparent  that  the  best 
results  depend  on  what  you  wish  to  accomplish.  Values  of  K3  =  2.0  and  t. 
=  0.4  seem  to  show  the  best  damped  control  responses,  but  values  of 
K3  =  1.0  and  tc  =  0.2  gave  the  best  accuracy.  The  remainder  of  the  LOS+ 
guidance  algorithm  tests  are  run  using  the  best  accuracy  results. 
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Figure  4-14  FOM  For  All  LOS+  Guidance  Runs 

4 . 7  LOS+  5  Degree  Off-Angle  Tests 

The  less  dynamic  5  degree  off-angle  shot  is  run  to  test  the  algorithm 
with  less  stringent  up-front  control  requirements,  and  to  test  the  effect 
of  a  5  G  aircraft  turn  on  the  guidance  situation.  The  bar  graph  labeled 
Sen  #1  and  Sen  #2  in  Figure  4-14  represent  these  two  tests.  Figure  4-15 
and  4-16  give  the  ground  map  and  control  graphs  of  these  two  runs.  In  the 
non-maneuvering  shot  of  Figure  4-15  some  erratic  behavior  in  the  heading 
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and  pitch  graph  appears  as  a  glitch  between  the  first  and  second  TOF 
points.  Although  the  data  in  the  range  is  confirmed  by  hand  calculations, 
an  error  or  a  cause  for  this  behavior  is  not  determined.  This  phenomenon 
shows  up  equally  in  the  azimuth  and  pitch  aspect  of  the  algorithm  leading 
suspicion  to  a  dynamic  pole  location  (as  seen  in  the  LOS  algorithm 
development)  as  the  cause.  This  is  supported  by  the  fact  that  this 
phenomena  changes  as  the  different  control  parameters  are  used  in  Figure 
4-10  through  4-13. 

In  Figure  4-16  one  can  see  the  effect  on  this  algorithm  when  the 
aircraft  maneuvers.  This  is  a  special  problem  for  this  algorithm  because 
it  causes  motion  in  the  future  LOS  to  the  target;  motion  that  is  not 
accounted  for  in  the  algorithm  development.  The  result  is  an  erratic 
flight  path  that  is  attempting  to  settle  on  an  unpredicted  future  line-of- 
sight .  Notice  also  the  pitch  channel  is  unchanged  from  the  previous  shot 
since  the  scenario  does  not  change  appreciably  in  pitch.  Although  this  is 
a  less  dynamic  control  challenge  the  results  of  this  algorithm  are 
unimpressive.  The  shots  are  useful  for  better  understanding  the  LOS+ 
guidance  algorithms  shortfalls. 

4 . 8  Radar  Off  LOS+  Tests 

When  the  sensor  model  uses  a  very  simple  3DOF  missile  model  to  estimate 
the  range  and  range  rate  measurements,  the  algorithm' s  dependence  on 
accurate  range  data  is  examined.  Figure  4-17  shows  the  no-radar  scenario 
and  it  should  be  compared  to  Figure  4-10  which  shows  the  same  shot  with 
accurate  range  and  range  rate  data.  Although  there  is  an  increase  in  the 
miss  distance  the  algorithms  performance  is  relatively  unaffected  by  poor 
range  and  range  rate  data.  Poor  range  data  causes  the  time  to  go 
calculation  to  suffer  but  if  the  guidance  commands  are  released  from  that 
dependence  the  results  can  be  relatively  independent  of  range  and  range 
rate.  (Much  the  same  as  the  LOS  guidance  was.)  Recall  that  the  time-to- 
go  accuracy  dependence  is  forced  to  zero  for  the  last  half  second  of 
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Figure  4-17  LOS+  Radar  Off  Test 


flight.  Although  this  algorithm  does  not  have  a  strong  range  and  range 
rate  dependence  it  has  enough  other  drawbacks  to  dwarf  this  advantage. 

4 . 9  Pursuit  Control  Parameter  Testing 

The  runs  used  to  establish  the  pursuit  algorithm  control  gain  are 
shown  in  Figure  4-18  through  4-22.  The  only  parameter  changed  for  each  of 
these  shots  is  Kj,  the  gain  for  the  pursuit  guidance  algorithm.  It  is 
typical  for  the  pursuit  algorithm  to  have  a  large  acceleration  towards  the 
end  of  its  flight.  This  is  caused  by  an  increase  in  the  missile-to- 
target  line-of-sight  rate  as  the  missile  range  to  the  target  decreases. 
In  Figure  4-18  one  can  see  this  high-G  turn  commanded  at  the  end  of  the 
missile  flight. 

An  interesting  aspect  of  this  “end  game"  turn  is  that  it  is  in  the 
opposite  direction  of  a  normal  pursuit  engagement.  In  the  normal  pursuit 
engagement,  i.e.  whare  a  missile  sensor  gives  the  LOS  to  the  target,  the 
missile  ends  up  pursuing  the  target,  and  flying  up  its  tail.  (Thus  the 
name  of  the  engagement . )  In  our  engagement  notice  that  the  tank  is  moving 
east  and  the  missile  requires  a  negative  heading  command  to  strike  the 
target;  bringing  it  into  the  nose  of  the  tank.  This  interesting  situation 
is  brought  about  by  the  fact  that  the  aircraft  is  doing  the  tracking,  not 
the  missile.  Thus,  in  relation  to  the  aircraft,  the  tank  is  moving  west 
and  the  pursuit  algorithm  is  performing  normally. 

In  Figure  4-19  through  4-22  are  some  runs  with  various  gains  to  show 
the  inability  to  correct  the  large  miss  distances  by  increasing  the  gain. 
Figure  4-23  shows  the  FOMs  for  all  the  pursuit  algorithm  shots.  Notice 
that  a  gain  of  K,  =  8  gives  the  best  performance  obtainable  with  this 
algorithm.  When  firing  a  hyper-velocity  missile  it  is  not  wise  to  put 
your  most  strenuous  turn  in  the  end  game  as  the  pursuit  algorithm  does. 
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Figure  4-1?'  Pursuit  Control  Parameter  Test  #2 
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Figure  4-21  Pursuit  Control  Parameter  Test  #4 
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4.10  Pursuit  5  Degree  Off-Anale  Test 

The  less  dynamic  5  degree  off-angle  shot  is  run  to  test  the  algorithm 
with  less  stringent  up-front  control  requirements,  and  to  test  the  effect 
of  a  5-G  aircraft  turn  on  the  guidance  situation.  The  bar  graph  labeled 
Sen  #1  and  Sen  #2  in  Figure  4-23  represent  these  two  tests.  Figure  4-24 
and  4-25  give  the  ground  map  and  control  graphs  of  these  two  runs. 

In  the  non  maneuvering  shot  of  Figure  4-24  one  can  still  see  a  large 
"end  game"  acceleration  commanded.  Again  this  is  caused  by  the  missile 
to  target  line  of  sight  rate  approaching  infinity  as  the  missile  closes  on 
the  target.  Remember  that  the  HVM  missile  is  a  hit  or  miss  "bullet”  and 
here,  even  with  no  sensor  errors  and  benign  launch  conditions,  the  missile 
only  passed  within  30  feet  of  the  target. 


4-31 


SCENARIO  GROUND  MAP 


4-33 


When  the  5-G  aircraft  maneuver  is  included  in  the  shot  the  miss 
distance  increases.  Comparing  Figure  4-25  with  4-24  gives  some  insight  to 
the  reason.  A  first  impulse  for  this  shot  is  that  it  is  independent  of 
an  aircraft  maneuver  since  the  aircraft  is  only  the  sensor  in  charge  of 
calculating  the  missile-to-target  line-of-sight  and  this  calculation 
should  be  relatively  independent  of  the  aircraft's  acceleration,  and  to 
some  extent  independent  of  its  position.  But  recalling  the  argument  that 
the  target  is  sensed  in  the  aircraft  reference  frame,  what  changes  here  is 
the  sensed  target  velocity.  Or,  in  reference  to  our  sensor  the  missile 
is  chasing  an  accelerating  target  and  this  causes  the  increased  miss 
distance.  One  can  see  the  increased  "end  game"  correction  attempt  when 
comparing  it  to  the  non-maneuvering  shot.  Other  than  this  perceived 
accelerating  target  problem  the  pursuit  algorithm  operates  the  same  for  a 
maneuvering  aircraft . 

4.11  Radar  Off  Pursuit  Tests 

Once  again  a  very  simple  3D0F  missile  model  is  used  to  estimate  the 
range  and  range  rate  measurements  to  examine  the  algorithm's  dependence  on 
accurate  range  data.  Figure  4-26  shows  the  no-radar  scenario  and  it 
should  be  compared  to  Figure  4-18  which  shows  the  same  shot  with  accurate 
range  and  range  rate  data.  Since  the  pursuit  algorithm  calculates  the 
missile  velocity  vector  direction  using  missile  range  and  range  rate,  i.e. 
see  equation  (3-21)  and  (3-22),  the  algorithm  performance  is  directly 
dependent  on  the  accuracy  of  this  data.  An  interesting  anomaly  can  be 
seen  in  Figure  4-26.  Notice  that  the  commanded  left  turn  increases  as 
the  missile  approaches  the  target.  Then,  the  range  estimator  says  the 
missile  has  past  the  target  and  the  command  abruptly  turns  the  missile 
into  a  right  turn.  This  last  right  turn  does,  coincidentally,  bring  the 
missile 
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towards  the  target,  but  not  because  the  algorithm  knows  the  missile  to 
target  line  of  sight.  It  has  been  unable  to  calculate  that  line  of  sight 
because  of  inaccurate  range  and  range  rate  data.  This  algorithm  needs  an 
accurate  range  and  range  rate  measurement  to  calculate  the  missiles 
heading  in  reference  to  the  target . 

4 . 12  Proportional  Control  Parameter  Testing 

The  various  runs  used  to  test  the  navigation  con.stant,  K3,  for  the 
proportional  guidance  algorithm  are  shown  in  Figure  4-27  through  4-31. 
The  only  guidance  parameter  changed  for  these  runs  is  Kj,  the  input 
variable  used  as  the  navigation  constant  when  the  proportional  algorithm 
is  selected.  As  before  in  these  figures  a  scenario  ground  map  that 

pictures  the  aircraft  track,  the  target  track  and  the  missile  flight  path 
is  included.  In  the  upper  left  corner  of  this  map  is  the  type  of 

algorithm  used,  some  figures  of  merit  and  the  control  parameters.  Inset 
into  each  of  the  scenario  maps  is  a  graph  of  the  heading  rate  commands  and 
graph  of  the  missile  heading  and  pitch. 

The  normal  navigation  constant  for  this  algorithm  ranges  from  2.5  to 
3.5  and  Figure  4-27  shows  test  results  with  K,  =  3.  This  shot  shows  a  very 

smooth  flight  path  that  takes  the  missile  to  an  interception  point  with 

the  target  as  designed.  The  miss  distance  of  2.7  feet  with  "perfect” 
single  precision  sensor  measurements  is  from  calculation  round-offs  and 
represents  an  equivalent  pointing  error  of  .26  mrad. 

Some  other  navigation  constant  runs  are  shown  in  Figures  4-28  through 
4-31  .  Notice  that  the  form  of  the  commanded  heading  rates  follows  a 
similar  pattern.  In  Figure  4-31  the  navigation  constant  of  Kj  =  6  is  used 
and  the  miss  distance  and  control  efficiency  FOMs  reach  another  minimum. 
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Figure  4-29  Proportional  Control  Parameter  Test  #3 
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Figure  4-31  Proportional  Control  Parameter  Test  #5 


Notice  here  that  more  of  the  turn  is  made  before  the  missile  reaches 
its  maximum  velocity.  The  higher  gain  also  puts  the  missile  on  a  good 
intercept  path  faster  than  the  lower  gain.  This  could  be  important  for 
shorter  range  shots  that  are  not  tested.  In  a  shorter  range  shot  the 
lower  gain  may  not  provide  the  correction  fast  enough  to  bring  the  missile 
in.  The  more  rounded  flight  path  of  Figure  4-27  is  selected  as  the 
nominal  control  parameter  because  it  is  the  more  conventional  navigation 
constant.  The  higher  gain  region  should  be  tested  in  shorter  range  in 
future  analysis. 

Figure  4-32  shows  all  the  FOM  results  for  the  proportional  algorithm 
runs.  The  first  five  graphs  show  the  tests  run  on  navigation  constant 
variations.  The  navigation  constant  value  of  3  shows  a  minimum  in  both 
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Figure  4-32  FOM  For  All  Proportional  Guidance  Runs 
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FOMs  and  the  remainder  of  the  proportional  guidance  algorithm  tests  are 
done  with  this  navigation  constant . 

4.13  Proportional  5  Degree  Off-Angle  Tests  The  less  dynamic  5  degree 
off-angle  shot  is  run  to  test  the  algorithm  with  less  stringent  up-front 
control  requirements,  and  to  test  the  effect  of  a  5-G  aircraft  turn  on  the 
guidance  situation.  The  bar  graph  labeled  Sen  #1  and  Sen  #2  in  Figure 
4-32  represent  these  two  tests.  Figure  4-33  and  4-34  give  the  ground  map 
and  control  graphs  of  these  two  runs.  Comparing  these  two  graphs,  the 
missile  to  target  line  of  sight  rate  is  relatively  independent  of 
accelerations  by  the  sensing  aircraft.  Recall  from  the  pursuit  algorithm 
that  the  missile  to  target  line  of  sight  is  dependent  on  aircraft 
acceleration,  where  it  appeared  to  the  algorithm  as  target  accelerations. 

Here  the  aircraft  acceleration  has  minimal  effect  on  the  FOMs  and  the 
precision  of  this  algorithm  rivals  the  LOS  algorithm. 

4.14  Radar  Off  Proportional  Tests 

When  the  sensor  uses  a  very  simple  3D0F  missile  model  to  estimate  the 
range  and  range  rate  measurements  the  algorithm's  dependence  on  accurate 
range  data  is  examined.  Figure  4-35  shows  the  no-radar  scenario  and  it 
should  be  compared  to  Figure  4-27  which  shows  the  same  shot  with  accurate 
range  and  range-rate  data.  It  is  obvious  from  examining  the  figure  that 
the  proportional  guidance  algorithm  has  a  direct  dependence  on  accurate 
range  and  range  rate  measurements.  This  dependence  is  rooted  in  the 
calculation  of  the  missile  to  target  line  of  sight  rate.  This  rate 
requires  an  accurate  assessment  of  the  missile  velocity  vector  which 
cannot  be  made  without  accurate  range  and  range  rate  data.  As  in  the 
pursuit  algorithm,  notice  the  abrupt  change  in  the  commands  when  the  range 
estimator  shows  the  missile  past  the  target.  Notice  also  that  despite  the 
poor  range  and  range  rate  data  the  initial  portion  of  the  flight  shows  the 
smooth  and  efficient  turn  into  an  intercept  approach. 
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Figure  4-33  Proportional  5  Degree  Off  Angle  Test 
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Figure  4-35  Proportional  Radar  Off  Test 


The  strong  dependence  on  range  and  range  rate  accuracy  may  be  the  drawback 
of  the  algorithm  for  the  HVM.  If  there  is  a  means  of  providing  this 
information  this  algorithm  has  great  potential  for  the  HVM  being  guided  by 
the  launch  aircraft. 

4.15  Algorithm  Comparisons 

The  performance  of  each  tested  algorithm  under  no-noise  conditions  is  the 
initial  indicator  of  its  worthiness  to  the  task  of  guiding  the  HVM  missile. 

The  algorithms  performance  will  only  degrades  under  realistic  .loise 
conditions.  The  less  robust  solutions  will  degrade  more  rapidly.  Table  4-1 
shows  the  FOMs  and  some  other  results  of  each  run  described  in  this  chapter. 

Comparative  bar  graphs  of  the  FOMs  are  shown  for  the  best  30  degree  off- 
angle  results/  in  Figure  4-36;  the  5  degree  off-angle  tests,  in  Figure  4-37; 
and  the  radar-off  shots,  in  Figure  4-38. 
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Figure  4-36  Comparison  of  All  30  Degree  Off  Angle  Shots 

By  examining  these  results  it  is  clear  that  the  LOS  guidance  algorithm  is 
superior  to  any  of  the  others  in  no-noise  performance.  A  close  second  in 
overall  performance  is  the  proportional  algorithm.  Its  only  drawback  is  the 
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Table  4-1  Figures  of  Merit  And  Parameters  For  Each  Run 


RUN 

TIME 

IMPCT 

(sec) 

MISS 

DIST 

...(ft) 

MISS 

X 

(ft) 

MISS 

Y 

(ft) 

MISS 

Z 

(ft) 

TERM 

VEL 

(fps) 

ENRGY  TC 

USED 

(lb  msec) (sec) 

K3 

(nu) 

LOS  1 

2.88 

4.83 

4.73 

0.25 

-0.95 

4908 

2.51 

0.100 

0.003 

LOS  2 

2.90 

4.18 

-2.8  9 

-3.02 

-0.05 

4842 

2.90 

0.080 

0.003 

LOS  3 

2.91 

5.18 

-2.70 

-4.41 

-0.39 

4809 

2 . 92 

0.100 

0 . 002 

LOS  4 

2.87 

1.57 

1.14 

-0.98 

-0.45 

4940 

2.27 

0.100 

0.004 

LOS  5 

2.84 

3.52 

2.89 

-1.76 

-0.99 

5026 

1.84 

0.200 

0.004 

LOS  S#1 

2.76 

1.72 

-0.15 

-1.62 

-0.56 

5277 

0.73 

0.100 

0.003 

LOS  S#2 

2.78 

4.37 

4.29 

-0.21 

-0.82 

5076 

0.99 

0.100 

0.003 

LOS  NoR 

2.86 

4.60 

4.50 

0.35 

-0.87 

4945 

2.23 

0.100 

0.003 

LOS+  1 

3.53 

26.94 

12.01 

-22.53 

-8.59 

2923 

6.54 

0.200 

1.000 

LOS+  2 

3.50 

53.04 

19.56 

-48.58 

-8.37 

3124 

9.73 

0.200 

2.000 

LOS+  3 

3.16 

106.65 

44.95 

-95.20 

-17.02 

3965 

5.78 

0.400 

4.000 

LOS+  4 

3.06 

107.99 

46.16 

-96.03 

-17.60 

4482 

4.41 

0.400 

2.000 

LOS+  S#1 

2.86 

25.98 

14.81 

-17.28 

-12.52 

4837 

1.43 

0.200 

1.000 

LOS+  S#2 

3. IF 

98.75 

77.25 

-60.80 

-9.31 

3753 

4.31 

0.200 

1.000 

LOS+  NoR 

3.46 

54.71 

21.29 

-49.46 

-9.69 

3074 

6.14 

0.200 

1.000 

PUR  1 

2.89 

104.57 

31.35 

-98.20 

-17.63 

4596 

3.71 

*  ★ 

4.000 

PUR  2 

2.87 

65.62 

18.32 

-62.20 

-10.04 

4647 

4.95 

★  ★ 

8.000 

PUR  3 

2.87 

44.84 

11.69 

-42.81 

-6.46 

4657 

8 . 90 

★  ★ 

16.000 

PUR  4 

2.87 

37.50 

14.04 

-34.28 

-5.84 

4640 

11.81 

★  ★ 

24.000 

PUR  5 

2.94 

34.34 

10.55 

-32.29 

-5.03 

4456 

15.01 

*  * 

32.000 

PUR  S#1 

2.76 

30.02 

17.52 

-18.85 

-15.46 

5096 

2.60 

*  * 

8.000 

PUR  S#2 

2.79 

111.53 

73.47 

-83.16 

-11.22 

4774 

4.98 

★  ■* 

8.000 

PUR  NoR 

2.89 

103.73 

-27.39 

95.78 

28.93 

4389 

36.72 

*  * 

8.000 

PRO  1 

3.35 

1.02 

-0.14 

1.00 

-0.14 

3526 

5.09 

if  if 

2.000 

PRO  2 

3.21 

2.68 

-2.59 

-0.61 

0.37 

3958 

3.33 

if  if 

3.000 

PRO  3 

3.12 

3.40 

3.16 

1.21 

-0.41 

4251 

3.07 

*  *■ 

4.000 

PRO  4 

3.00 

2.58 

2.34 

1.05 

-0.27 

4637 

2.45 

*  * 

6.000 

PRO  S#1 

2.82 

1.03 

0.61 

0.82 

-0.14 

4972 

1.03 

*  * 

3.000 

PRO  S#2 

2.82 

1.26 

0.79 

0.97 

-0.18 

4972 

1.06 

★  it 

3.000 

PRO  NoR 

2.89 

>500.00 

*  * 

*  * 

*  * 

4389 

9.41 

it  ★ 

3.000 

strong  dependence  on  accurate  missile  range  and  range-rate  measurements. 
From  Figure  4-37  one  can  see  that  it  is  more  resilient  to  aircraft 
acceleration  effects  than  is  the  LOS  algorithm. 

The  effort  to  project  an  aircraft  position  in  time  and  fly  the  missile  to 
a  future  line  of  sight  is  a  difficult,  dynamic  control  problem  that  will 
require  a  varying  gain  control  for  effective  use.  The  LOS+  algorithm  with  a 
constant  gain  represents  a  control  nightmare  that  provides  neither  accuracy 
nor  efficiency. 

The  pursuit  algorithm  is  unsuitable  for  any  high  closing  speed  missile. 
Its  inefficiency  and  lacking  accuracy  for  the  HVM  is  prominent  in  the  graphs. 
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Figure  4-37  Comparison  of  All  5  Degree  Off  Angle  Shots 
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Figure  4-38  Comparison  of  All  Radar  Off  Shots 


V.  Conclusions  and  Recommendations 


5 . 1  Introduction . 

After  working  closely  with  the  proposed  HVM  weapon  system  for  the 
initial  testing,  an  insight  into  the  problem  and  the  "personality"  of  the 
algorithms  is  developed.  This  insight  provides  the  basis  for  more 
detailed  examination  of  selected  promising  algorithms  before  the  HVM  can 
be  effectively  steered  by  the  launch  aircraft.  Although  this  work  has 
been  in  a  no-noise  environment  a  good  understanding  of  the  physical 
limitations  of  each  algorithm  can  be  very  useful  in  selection  of  the  HVM' s 
guidance  scheme.  The  robustness  of  each  algorithm  can  only  be  tested  with 
a  sensitivity  analysis  that  includes  sensor  and  system  noise.  This 
chapter  states  the  conclusions  of  the  completed  effort  and  provides 
recommendations  for  future  study. 

5 . 2  Conclusions . 

Currently  used  missile  guidance-concepts  are  applicable  to  the  task  of 
guiding  a  hyper-velocity  missile  without  a  target  sensor  into  a  direct 
collision  course  with  a  ground  target.  The  following  conclusions  are 
drawn  from  the  no-noise  tests  accomplished  on  the  four  algorithms.  In 
each  case  it  must  be  considered  that  algorithm  robustness  and  stability  is 
not  adequately  tested  with  this  effort,  only  the  physical  constraints  of 
the  algorithms  within  the  HVM  weapon  system  geometry. 

5.2.1  LOS  Guidance  Algorithm.  When  the  LOS  guidance  algorithm  is 
developed  and  controlled  with  properly  tuned  control  parameters,  it  can 
effectively  guide  the  missile  into  the  target  under  a  variety  of  launch 
conditions .  The  LOS  algorithm  can  perform  well  with  only  modeled  range 
and  range-rate  data  because  it  is  relatively  insensitive  to  these 
parameters.  Thus  a  FLIR  sensor  providing  angle  data  on  the  missile  and 
target  is  adequate  for  this  algorithm.  Although  the  LOS  algorithm 
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represents  a  time  dependent,  non-linear  control  system,  a  constant  gain 
setting  can  provide  stable  and  accurate  results  under  the  no-noise 
simulation  accomplished  in  this  effort.  This  algorithm  preforms  the  best 
overall,  but  it  is  a  dynamic  control  system  that  must  be  tuned  in  the 
presence  of  system  noise  before  a  final  conclusion  can  be  drawn. 

5.2.2  LOS+  Guidance  Algorithm.  The  LOS+  algorithm,  devised  to  increase 
LOS  algorithm  efficiency  and  to  move  the  missile  plume  away  from  the 
aircraft  to  target  line-of-sight,  requires  time  varying  gain  control  to 
work  effectively.  Although  it  does  move  the  missile  plume  off  the  target 
line-of-sight  effectively,  the  non-optimal  control  and  its  dependence  on 
predicting  the  future  are  insurmountable  problems  in  this  effort.  A  gain 
balance  that  optimally  moves  the  missile  on  to  a  future,  stable  line-of- 
sight  and  provides  the  end  game  gain  necessary  for  accuracy  cannot  be 
obtained  with  a  constant  gain  as  attempted  with  this  effort.  If  this 
technique  will  work  it  will  require  a  time  variable  gain,  preferably  a 
optimal  guidance  solution. 

5.2.3  Pursuit  Guidance  Algorithm.  The  pursuit  algorithm  cannot  perform 
well  in  a  high  closing  speed  environment.  The  HVM  presents  an 
exceptionally  high  closing  speed  environment  and  the  pursuit  algorithm 
calls  for  end  game  accelerations  that  cannot  be  met  by  the  missile. 
Although  this  guidance  principle  is  very  popular  with  air-to-air 
intercepts,  and  with  some  missile  IR  seekers  where  an  up  the  tail  shot  is 
essential,  it  is  not  a  candidate  for  the  HVM  environment. 

5.2.4  Proportional  Guidance  Algorithm.  The  proportional  guidance 
algorithm  performs  admirably  in  the  HVM  weapon  system.  It  shows  excellent 
independence  of  launch  aircraft  maneuvers  and  with  accurate  sensor  data  it 
can  deliver  the  accuracy  required  for  a  direct  target  hit.  The 
proportional  guidance  algorithm  is  more  stable  than  the  LOS  algorithm  and 
may  out  perform  this  algorithm  in  the  presence  of  system  noises.  The 
proportional  algorithm  has  a  strong  dependence  on  aircraft  to  missile 
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range  and  range-rate  measurement  accuracies.  The  utility  of  this 
algorithm  will  depend  on  the  AFTI/F-16's  ability  to  provide  accurate  range 
and  range-rate  information. 

5.2.5  Conclusion  Summary.  This  study  effort  is  the  initial  examination 
of  guidance  concepts  for  the  HVM  missile.  The  no-noise  computer 
simulation  runs  only  establish  a  baseline  performance  of  the  algorithms. 
Additional  work  is  required  to  determine  the  algorithm's  performance  when 
sensor  inaccuracies  and  measurement  noise  are  added  to  the  simulation. 
From  the  noiseless  results  the  algorithms  which  provide  acceptable 
performance  are  the  I.OS  guidance  algorithm  and  the  proportional  guidance 
algorithm. 

5 . 3  Recommendations 

There  are  four  areas  where  recommendations  are  made  with  confidence. 
Based  on  the  results  of  this  work  and  the  author's  familiarity  with  the 
algorithms  the  best  performing  algorithms  are  selected.  The  effect  of  the 
missile  plume  interference  with  the  FLIR  requires  further  examination,  but 
an  approach  to  solve  this  potential  problem  is  included.  Lastly,  some 
recommendations  for  the  next  logical  effort,  a  sensor  error  analysis,  are 
included.  Each  of  these  recommendations  is  discussed  in  the  following 
paragraphs . 

5.3.1  Algorithm  Selection.  If  the  HVM  development  continues  with  a 
missile  that  can  only  sense  rotation  and  execute  azimuth  and  pitch  turn 
commands,  either  the  LOS  or  the  proportional  guidance  algorithm  should  be 
developed  for  it .  If  studies  indicate  that  a  radar  can  not  provide 
accurate  range  and  range-rate  data  for  the  HVM,  or  if  tactical  users 
continue  to  expr,-3S  a  strong  desire  to  launch  with  out  emitting  (i.e. 
radar  off) ,  then  the  LOS  algorithm  is  best  suited  to  guide  an  HVM  to  a 
direct  hit  on  a  ground  target. 

5.3.2  Optimal  Guidance  Algorithm  Development.  Several  studies  indicate 
that  the  proportional  guidance  algorithm  preforms  as  well  as  an  optimal 
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guidance  solution  in  most  situations.  [6:237-238]  In  the  tests  shown  here 
with  no  system  noise,  the  proportional  guidance  algorithm  does  not 
significantly  out  perform  the  LOS  algorithm.  If  an  optimal  guidance 
scheme  is  pursued  for  the  HVM  it  should  be  centered  on  the  final  line-of- 
sight  concept  of  the  LOS+  algorithm.  This  technique  can  still  provide  the 
high  accuracy  and  range  independence  of  the  LOS  algorithm,  while  moving 
the  missile  off  the  launch-aircraft  to  target  line-of-sight . 

5.3.3  Missile  Plume  Obstruction  Problems.  If  it  is  determined  that  the 
missile  plume  obstructs  the  target  tracking,  or  that  the  target  FLIP, 
signature  hinders  tracking  the  coasting  missile,  a  hybrid  LOS/Proportional 
guidance  algorithm  should  be  evaluated.  The  proportional  algorithm  would 
take  the  missile  off  the  line-of-sight  during  rocket  burn,  then  the  LOS 
algorithm  can  give  the  final  accuracy  and  range  independence  necessary  for 
the  direct  hit.  This  can  be  done  by  calculating  the  commands  for  both 
algorithms  and  time  weight  them  before  sending  a  command  set  to  the 
missile.  This  would  represent  a  simpler  solution  than  a  optimal  guidance 
algorithm  development,  but  would  need  to  be  validated  in  a  simulation  and 
compared  to  the  other  algorithms. 

This  hybrid  solution  combines  two  nearly  optimal  solutions,  gives  a 
controllable  time  where  the  missile  is  off  the  subject  line-of-sight,  and 
retains  the  robust  accuracy  of  the  LOS  algorithm.  Examining  the  30  degree 
off-angle  shots  for  both  proportional  and  LOS  guidance  shows  that  the 
proportional  guidance  lines  up  the  missile  for  a  final  LOS  guidance 
scheme.  In  fact  each  of  the  shots  tested  shows  this  hybrid  to  be  a 
feasible  solution. 

5.3.4  Suggestions  For  Further  Study.  Several  areas  are  suggested  for 
further  study  of  these  algorithms  in  the  no-noise  situation  in  which  they 
were  tested.  The  proportional  algorithm  shows  a  range  dependence  in  this 
test  but  a  sensitivity  to  range  and  range-rate  biases  can  be  further 
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explored  with  this  simulation  capability.  Sensitivity  to  FLIR  angle 
biases  can  also  be  explored  with  this  no-noise  simulation. 

The  miss  distances  used  for  a  FOM  in  this  study  are  most  likely  caused 
by  data  latency.  In  the  time  between  each  missile  update  at  50  Hz,  the 
missile  moves  about  100  feet  closer  to  the  target.  The  resulting  piase 
lag  likely  causes  the  terminal  oscillations  and  larger  miss  distances. 
The  accuracy  of  range  and  angle  measurements  provided  to  the  guidance 
algorithm  is  limited  in  this  simulation  by  the  24  bit  mantissa  word 
length.  So  the  miss  distances  are  a  function  of  the  guidance  algorithms 
ability  to  control  the  hyper-velocity  projectile,  the  latency  in  the  data 
updates  and  the  missile's  second  order  response  to  a  command.  A  command 
rate  study  can  evaluate  the  missile  command  update  time  requirements  and 
the  effects  of  missed  updates  on  the  missile  accuracy.  The  performance 
sensitivity  to  sensor  update  times,  or  filter  update  times  with  different 
missile  update  times  can  also  be  explored  with  this  model.  This  study  can 
assist  the  design  of  the  update  scheme  used  on  the  HVM  missile,  its 
robustness,  and  its  timing  requirements. 

The  LOS  algorithm  was  analyzed  using  small  angle  assumptions  to  make 
alinear  control  block  diagram.  A  functional  analysis  of  this  non-linear 
algorithm  can  provide  insight  to  constructing  a  time  varying  gain  control 
for  this  approach.  The  3DOF  no-noise  model  constructed  for  this  effort 
can  be  used  in  several  of  these  study  areas. 

5.3.4  Sensor  Error  Analysis.  As  tie  HVM  design  matures  and  emphasis  is 
placed  on  its  development,  the  results  of  this  effort  should  be  expanded 
to  include  a  sensor  error  analysis.  The  ability  of  a  FLIR  to  track  a  HVM 
missile  throughout  the  course  of  its  flight  needs  to  be  analyzed.  Capp' s 
Kalman  filter  development  [4:55]  and  the  MSOFE  Kalman  filter  modeling 
tools  should  then  be  used  to  perform  a  system  simulation  with  realistic 
noise.  The  tests  should  continue  with  the  LOS  algorithm  and  perhaps  the 
LOS/Proportional  hybrid  described  in  paragraph  5.3.3  above.  The  results 
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of  such  a  test  are  expected  to  continue  showing  the  superiority  of  the  LOS 
guidance  algorithm  for  the  HVM  weapon  system. 
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Various  Additional  Scenerio  Launches 
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The  following  additional  sceneries  were  explored  to  test  the  conclusions 
of  this  effort.  The  data  is  presented  here  for  completeness  and 
additional  insight  into  selected  guidance  algorithms. 

Table  A-1  Missile  Parameters  Used 


Description 

Default 

Value 

Unit 

Coded 

Variable 

- 

MSL  THRUST  TIME 

1.5 

SEC 

TBURN 

2 

MSL  THRUST 

5000 

LBS 

THRUST 

MSL  BODY  WEIGHT 

25.9 

LBS 

WTBD 

4 

MSL  FUEL  WEIGHT 

32.17 

LBS 

WTFU 

MSL  DIAMETER 

3.8 

IN 

MDIA 

MSL  LENGTH 

7.0 

FT 

MLNGTH 

"7 

MSL  Wn 

44.0 

RPS 

WN 

P 

MSL  DAMPING 

.7071 

NU 

ZETA 

9 

MSL  MAX  G 

102.0 

G 

GMAX 

10 

COEFF  LIFT  DRAG 

1 .74E-4 

NU 

KDRAG 

1  1 

MAX  MSL  FLIGHT 

TIME  4.0 

SEC 

TMAX 

Table  A-2  Scenario  Parameters  Used 

Default 

Coded 

Description 

Value 

Unit 

Variable 

AIRCRAFT  VEL 

650.00 

FPS 

VAC 

2 

ACFT  HEADING 

60.0 

DEG 

HDGA 

ACFT  ALTITUDE 

1000.0 

FT 

AALT 

4 

ACFT  LAT  ACCEL 

000.0 

G 

AACf 

C 

ACFT  ACCEL  PITCH 

00.0 

DEG 

APITCH 

6 

RANGE  TO  TARGET 

5280.0 

FT 

RT 

7 

LOS  TO  TARGET 

30.0 

DEG 

LOST 

8 

TARGET  VEL 

50.0 

MPH 

VTGT 

9 

TARGET  HEADING 

90.0 

DEG 

HDGT 

10 

TARGET  LAT  ACCEL 

0.0 

FPS^ 

2  ATGT 

X  a. 

SYSTEM  INTEG  TIME 

.002 

SEC 

DT 

12 

MSL  UPDATE  TIMES 

.02 

SEC 

DTMU 

13 

SENSOR  UPDATE  TIMES 

.02 

SEC 

DTSU 

14 

OUTPUT  TIMES 

.02 

SEC 

DTOUT 

1  Z 

FILTER  UPDATE  TIME 

.02 

SEC 

DTKU 
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Figure  A— 3  One  Mile  Shot 


SCENARIO  GROUND  MAP 


EAST  In  Thousand  Feet 


<3  =  Q.Q03 
TC  =  0.10 

1  Mile  Shot,  No  Radar,  Acft 
Tgt  Hdg  -90  Deg,  +5G  Turn 


o» 


O 


SCENARIO  GROUND  MAP 


(0 


A-8 


EAST  In  Thousand  Fee 


o 

J3 

U 

0) 

rH 

■H 

X 


01 

C 

O 


(V 

u 

3 

O' 

•H 


A- 9 
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Table  A-3  Missile  Parameters  Used 


Description 

Default 

Value 

Unit 

Coded 

Variable 

1 

MSL  THRUST  TIME 

1.5 

SEC 

TBURN 

2 

MSL  THRUST 

5000 

LBS 

THRUST 

MSL  BODY  WEIGHT 

25.9 

LBS 

WTBD 

4 

MSL  FUEL  WEIGHT 

32 . 17 

LBS 

WTFU 

c; 

MSL  DIAMETER 

3.8 

IN 

MDIA 

6 

MSL  LENGTH 

7.0 

FT 

MLNGTH 

7 

MSL  Wn 

44.0 

RPS 

WN 

8 

MSL  DAMPING 

.7071 

NU 

ZETA 

9 

MSL  MAX  G 

102.0 

G 

GMAX 

10 

COEFF  LIFT  DRAG 

1 .74E-4 

NU 

KDRAG 

11 

MAX  MSL  FLIGHT  TIME  4.0 

SEC 

TMAX 

Table  A-4  Scenario  Parameters  Used 

Default  Coded 

Description _ Value _ Unit  Variable 


1 

AIRCRAFT  VEL 

650.00 

FPS 

VAC 

O 

ACFT  HEADING 

90.0 

DEG 

HDGA 

ACFT  ALTITUDE 

1000.0 

FT 

AALT 

4 

ACFT  LAT  ACCEL 

000.0 

G 

AACf 

5 

ACFT  ACCEL  PITCH 

00.0 

DEG 

APITCH 

6 

RANGE  TO  TARGET 

5280.0 

FT 

RT 

LOS  TO  TARGET 

30.0 

DEG 

LOST 

8 

TARGET  VEL 

50.0 

MPH 

VTGT 

9 

TARGET  HEADING 

-90.0 

DEG 

HDGT 

10 

TARGET  LAT  ACCEL 

160.0 

FPS''2 

ATGT 

11 

SYSTEM  INTEG  TIME 

.002 

SEC 

DT 

12 

MSL  UPDATE  TIMES 

.02 

SEC 

DTMU 

13 

SENSOR  UPDATE  TIMES 

.02 

SEC 

DTSU 

14 

OUTPUT  TIMES 

.02 

SEC 

DTOUT 

15 

FILTER  UPDATE  TIME 

.02 

SEC 

DTKU 
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